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EXECUTIVE SUMMARY 

Despite the outcome of Operation Desert Shield/Storm, post-operation analysis 
revealed a shortcoming in USAF airlift capability. The analysis showed that early 
presence was not sustainable and an early Iraqi attack could have inflicted more coalition 
casualties. This finding and the desire to gain insights into the mobility system prompted 
Congress to sponsor the Mobility Requirement Study (MRS). 

This thesis involves an optimization modelling of strategic airlift assets (aircraft 
and airfields) and is sponsored by the USAF Studies and Analyses Agency. It is an 
outgrowth of the USAF and Joint Staff work for MRS. The thesis attempts to provide 
insights into the strategic airlift system. Basically, it attempts to provide answers to 
airlift mobility questions such as, what are the system "bottle-necks" and which set of 
alternative strategic airlift fleets is the most effective for the scenario. 

In this research, the strategic airlift assets optimization problem is formulated as 
a multi-commodity, multi-period network flow model with a large number of side 
constraints using linear programming (LP). The model, implemented on the General 
Algebraic Modelling System (GAMS), minimizes late deliveries subject to system 
constraints such as aircraft utilization rate and aircraft handling capacity of an airfield. 

This optimization model can provide broad analytic insights into the strategic 
airlift system. For any particular set of input data, the insights gained include answers 
to the following mobility questions: 1) Are the assets adequate for the scenario? 2) 


What are the impacts of a shortfall in airlift capability? 3) What are the system "bottlenecks"? 
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This thesis demonstrates that such an LP model has sufficiently fast response time 
to be a viable planning tool in today's political environment, where major regional 
conflicts can emerge quickly and simultaneously. The model can lend support to the 
study of military options at the planning and acquisition stages, as well as enable planners 


to quickly assess the impact of any shortfall in airlift capability. 
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I. INTRODUCTION 
In October 1993, the Naval Postgraduate School (NPS) Operations Research 
Department was formally enlisted by the United States Air Force Studies and Analyses 
Agency (USAF/SAA) to develop a Strategic Airlift Assets optimization model. This 
model would be used as a planning tool in analyzing mobility questions such as the 


impact of any shortfall in airlift capability for a particular major regional contingency. 


A. BACKGROUND 

Despite the outcome of Operation Desert Shield/Storm, post-operation analysis 
revealed a shortcoming in USAF airlift capability. The analysis showed that early 
presence was not sustainable and an early Iraqi attack could have inflicted more coalition 
casualties. It is for this reason that Congress commissioned a Mobility Requirement 
Study (MRS) in 1991. This comprehensive study examined all aspects of the mobility 
question, from domestic transportation. to both inter-theater and  intra-theater 
requirements. 

Two linear program (LP) optimization models arose directly from the 
Congressionally-mandated study and form the background of this thesis. They are the 
Joint Staff's Force Structure Resource, and Assessment Directorate (J8)'s Mobility 
Optimization Model [Ref. 1] and the USAF Studies and Analyses Agency's Thruput 


Model [Ref. 2]. 


l. Mobility Optimization Model 

The Mobility Optimization Model (MOM), a product of MRS, was developed 
by J8 with assistance from NPS within a short period of six weeks. MOM was 
developed to aid in determining the proper level and mix of lift assets (air and sea) 
necessary to support US power projection needs into the 21st century. Although MOM 
serves its purpose well, it is not suitable for answering the mobility questions currently 
sought by USAF/SAA because it does not model airlift constraints in sufficient detail; 
e.g., the maximum number of planes on the ground (MOG) constraint at intermediate 
airfields is not directly modelled. Some of MOM's technology has, however, been 


adopted by this thesis. 


2. Thruput Model 

The static Thruput optimization model was developed by USAF/SAA to help 
evaluate alternative military strategic airlift fleets, route structures and basing schemes 
for different scenarios. Although the Thruput model has proven to be useful, it has 
inherent limitations due to its static nature (single period). For example, important 
constraints such as delivery time windows cannot be modelled when the time domain is 
not incorporated. In addition, the model is unable to provide answers to key questions 
that concern decision makers such as: On what day will unit x commence movement and 


on what day will it be closed? 


B. THESIS GOAL 

The goal of this thesis is to develop a multi-period Strategic Airlift Assets LP 
optimization model using the General Algebraic Modelling System (GAMS) [Ref. 3] for 
the USAF. This new optimization model is an outgrowth of the USAF and Joint Staff 
work for MRS. It attempts to combine the best features of both models by retaining 
Thruput’s level of system detail while incorporating MOM’s treatment of the time 
domain. Specifically, the new model strives to provide answers to mobility questions, 


some of which are not available from existing models. For example: 


@ Are the given sets of aircraft and airfields adequate for the deployment scenario? 


@ If the assets are inadequate, what is the impact of the shortfall in terms of tons of 
equipment and number of men delivered to the theater late? 


e Which of the assets (aircraft or airfield) are "bottle-necks" in the system? Is there 
a need to increase the number of aircraft (reserves or civilian aircraft) or a need 
to negotiate for access to more foreign airfields? 


e Which set of alternative strategic airlift fleets is the most effective for the scenario? 


e On what day will unit x commence movement and on what day will it be closed? 


C. PROBLEM STATEMENT 

The precise nature of the airlift problem is: Given the movement requirements for 
both troops and equipment for a particular scenario, and the aircraft and airfield assets 
available, find the optimal combination of airlift mission assignments by number and type 
of aircraft for each unit, routing structure, airlift mission start time and cargo type to 


carry. The objective of the model is to minimize penalties (weighted by each unit's 


priority) for late deliveries and undelivered cargo, subject to system constraints such as 


bounds on aircraft utilization rates and the aircraft handling capacity of each airfield. 


D. JUSTIFICATION FOR MODEL 
The justifications for developing the Strategic Airlift Assets optimization model are 


the need to optimize and the need for quick answers. 


€ Need to optimize: With shrinking military budgets, the number of both aircraft 
and airfields available are set to decline. Further, the desire to have early force 
closure so as to reduce the vulnerability period when US forces may be 
outnumbered imposes great strains on the mobility system. These two factors 
imply that inevitably, demand will exceed supply and there is a need to optimally 
utilize both aircraft and airfield assets. 


€ Need to have quick answers: Today's political environment has great uncertainties 
in which a major regional conflict can emerge very quickly (for example Somalia). 
In this respect, an optimization model such as the strategic airlift model which has 
good turnaround time would be a viable and useful planning tool.! In fact, its 
quick turnaround time may allow the optimization model to be used in a 
complementary manner with other models which have different virtues; e.g., a 
simulation model can capture more detail but generally takes longer to run. One 
possible complementary mode of operation is to use the outputs of the optimization 
model as inputs for the simulation model, e.g., number of aircraft and airfield 
assets to use. 


'The current optimization model has a total response time of under 2 hrs for a typical 
size problem. It takes approximately an hour and a half to enter new data and the LP 
is generated and solved in under 2 minutes on an IBM RS6000 model 590 workstation. 


II. OVERVIEW OF MODEL 
This chapter gives a broad overview of pertinent features of the airlift system, 
modelling assumptions, and a verbal description ("word formulation") of the objective 
function and constraints. Also included in this chapter is a brief discussion on the origin 


of data sources, the need for data aggregation and model reduction techniques. 


A. MODEL FEATURES 

The airlift system has its own peculiar features and modes of operation. Most of 
these have been incorporated in the model to make it as representative of real time 
operation as possible. Others, such as the use of tanker aircraft for aerial refueling of 
airlift aircraft are recommended as follow-on work. One possible way to incorporate the 
use of tankers is the "airbase in the sky" concept utilized by RAND’s CONOP 
optimization model [Ref. 4]. The major aspects of the airlift system captured by the 


model include: 


@® Multiple origin/destination airfields (source and sink nodes): This feature is 
representative of the airlift system that typically utilizes multiple origin, enroute 
and destination airfields. The model can facilitate the study of simultaneous major 
regional contingency plans. 


e Flexible routing structure: The air route structure supported by the model includes 
delivery/recovery routes with zero to three enroute stops. This provision allows 
for use of a mixture of both short-range and long-range aircraft. It also allows the 
LP to choose between higher-payload, shorter-range flights (more enroute stops and 
more ground time) and lower-payload, longer-range flights. In addition to the 
variable number of enroute stops, the model also allows the same aircraft to fly 
different delivery and recovery routes for improved realism and efficiency. 


€ Aircraft-to-route restriction: Provisions have been made to allow the user to 
impose aircraft-to-route restrictions; e.g., military aircraft may use military 
airfields for enroute stops. This particular provision arises because Air Mobility 
Command (AMC) may rely on civilian airliners to augment USAF aircraft in a 
deployment under the Civil Reserve Airfleet (CRAF) program. It is necessary to 
distinguish between the USAF and CRAF aircraft because of possible route 
restrictions. For example, military aircraft may not be permitted to land in certain 
civilian airfields. The model also allows the user to route aircraft through specific 
recovery bases, as in AMC'S typical concept of operations. 


€ Aircraft assets can be added over time: This adds realism to the model since 
CRAF and reserves aircraft typically take time to mobilize and are generally not 
available until later in the deployment. 


e Delivery time windows: In a deployment, a unit is ready to move on a certain date 
and has to arrive at the theater by a specific date; 1.e., there is a Time Phase Force 
Deployment Date (TPFDD) associated with each unit. This aspect of the problem 
has been incorporated in the model through user-specified dates for each unit. 


B. ASSUMPTIONS 
Listed below are the major assumptions made by the model. These assumptions 


are made due to the nature of the data available or to avoid computational intractability. 


e Inventoried aircraft at origin/destination airfields do not affect the aircraft handling 
capacity of the airfield: This assumption is made in the modelling of the aircraft 
handling capacity constraint of each airfield. It is not strictly valid since an 
inventoried aircraft takes up parking space even if it is not consuming services. 
However, since the aircraft handling capacity figures are based on a host of other 
factors like material handling equipment (MHE), servicing and fuel availability, it 
is probably more accurate to model the problem as such. Segregation of airfield 
Capacity data (into parking spaces and services) is necessary before the constraints 
can be modelled appropriately. This recommendation has been made to the USAF. 


@ Deterministic ground time: Aircraft turnaround times for onload/offload of cargo 
and enroute refuelling are assumed to be known constants, although they are 
naturally stochastic. This assumption will result in an optimistic solution from the 
LP as time deviations could cause some aircraft to wait on the ground while others 
are being serviced. The stochastic nature of ground time is however, offset 


somewhat by a MOG efficiency factor Introduced in the formulation to soften the 
Impact of randomness. Further research and validation are required to assess the 
Impact of this assumption and the offsetting technique. 

e No restriction on airfield operating hours: The assumption that all airfields are 
twenty-four hour capable is included as a caveat for future modelers because an 
airfield may have operating hours restriction. Airfield operating hours do not 
affect the current model since the time resolution (in days) does not warrant a 
discrimination of airfield operating hours as arrivals at each airfield is rounded to 
the nearest day; i.e., no distinction between day and night arrival etc. Airfield 
operating hours may become important if the time resolution of the model is 
reduced to say six hour blocks. 


C. "WORD FORMULATION" 
This section gives a verbal description of the objective function and constraints of 
the Strategic Airlift Assets optimization model. Discussions on the mathematical 


formulation and implementation details can be found in Chapters III and IV respectively. 


1. Objective Function 
The purpose of the optimization model is to determine a schedule of airlift 
missions with troops and cargo to carry, and recovery missions, so as to minimize the 
total weighted penalties of late deliveries and undelivered cargo (weighted by the 


movement priority of a unit). 


2. Constraints 
The constraints of the airlift system modelled can be broadly grouped into the 
following categories: Demand satisfaction, aircraft balance, aircraft physical limitation, 


aircraft utilization rate and aircraft handling capacity of an airfield. 


e Demand Satisfaction Constraints: Demands on the airlift system are the movement 
requirements for troops and cargo. The demand satisfaction constraints for troops 
try to ensure that the number of troops that arrive within the permitted time 
window for each unit is greater than the number required, with shortfalls accounted 
for by deviation variables that are penalized in the objective function. The demand 
satisfaction constraints for each unit's cargo serve the same purpose for the 
required equipment. 


€ Aircraft Balance Constraints: These constraints keep physical count of aircraft by 
type in each time period. They ensure that the aircraft assets are used only if and 
when they are available. 

€ Aircraft Physical Limitation Constraints: There are three different kinds of 
constraints on the physical limitations of aircraft. These constraints ensure that 
limits on troop carriage capacity, maximum payload, and cabin floor space 
constraints are observed at all times. 

€ Aircraft Utilization Rate Constraints: These constraints ensure that the average 
flying hours consumed per aircraft per day is less than AMC's established 
utilization rate for each aircraft type. 

e Aircraft Handling Capacity of Airfield (MOG): These constraints ensure that the 


number of aircraft routed through an airfield each day can be handled by the 
airfield. 


D. DATA SOURCES, AGGREGATION AND MODEL REDUCTION 


1. Data Sources 
The data inputs used in trial runs of the model to date were obtained through 
the assistance of USAF/SAA. Most of these data are scenario-specific and can be found 
in documents and planning guides. 
For example, the airfields’ aircraft handling capacity figures (or MOG 


figures) were developed by AMC for the Mobility Requirements Study and data such as 


space or maximum payload of an aircraft for a particular range were obtained from 


aircraft technical manuals. Appendix A gives a listing of the data sources. 


2. Data Aggregation 

The need for data aggregation as well as efficient handling of multi- 
dimensional modeling entities was envisaged early during model development. For 
example, there is a set of decision variables with eight indices that would have over 9 
millions members if all combinations are included.” This set must be reduced drastically 
in model implementation in order for the LP to be tractable. For this reason, the 
geographic and requirements data were reviewed thoroughly with USAF/SAA and 
aggregated whenever reasonable to do so. For example, airfields in close proximity (up 
to about 300 nautical miles) are collapsed into a single entity with their MOG figures 
aggregated, and units with similar origins, destinations and time windows are represented 


as aggregate movement requirements. 


3. | Model Reduction 

In addition to data aggregation, logical conditional checks are used 
extensively in the model to reduce the number of decision variables. The conditional 
checks serve to eradicate those combination of decision variables that are known to have 
zero values at optimality and hence help to reduce the problem size. For a simple 
example, if there are no Boeing 747 aircraft (CRAF) available until day 10, then all 
“The variable has one index each for unit, time period and aircraft type, and 5 
indices representing airfields along a route. For a problem with 30 time periods, 20 


units, 5 aircraft types and 5 members per airfield index, with all possible index 
combinations counted, there would be 9,375,000 variables. 


airlift missions associated with the B747 fleet for the first 9 days are removed from the 
LP formulation. Unnecessary constraints are eliminated in a similar manner. 

These logical checks and reductions are implemented in a generic fashion, so 
in contrast to the aggregations they do not have to be reconsidered when the user changes 


the data. 
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IM. OPTIMIZATION MODEL 
This chapter gives a mathematical formulation of the conceptual optimization model 


discussed in Chapter II. Implementation aspects of the model are covered in Chapter IV. 


A. METHODOLOGY 

The strategic airlift assets optimization problem is basically formulated as a multi- 
commodity, multi-period network flow model with a large number of side constraints. 
Two key concepts are employed in the model. The first key concept 1s the use of a time 
index to track the locations of aircraft for each time period. Knowing when an aircraft 
will arrive at a particular airfield allows the aircraft handling capacity constraint of an 
airfield to be modelled appropriately. The second key concept is model reduction 
through data aggregation and the removal of unnecessary decision variables and 


constraints prior to LP optimization. 


B. MATHEMATICAL FORMULATION 


1. Indices 


unit; e.g., 82nd Airborne 
aircraft type; e.g., C5, C141 
t time period in days 
generic airfield; 1.e., origin, enroute and destination airfields 
origin airfield 
destination airfield 
route 


EE ca n UM 


11 


2 Index Sets 


a. Airfield Index Sets 
set of available airfields 


ISB origin airfields 
KSB destination airfields 


b. Aircraft Index Sets 


Α set of avallable aircraft types 

Ay S A aircraft capable of hauling bulk-sized cargo 
Aver S Apu aircraft capable of hauling over-sized cargo 
Aou S Ao aircraft capable of hauling out-sized cargo 


C. Route Index Sets 


R set of available routes 

RER permissible routes for aircraft a 

R,, SR permissible routes for aircraft a that utilizes airfield b 

R, SR permissible routes for aircraft a that connect the (i,k) origin- 
destination pair 

DR. R delivery routes that originate from origin i 

RR, €R recovery routes that originate from destination K 


d. Time Index Sets 
T set of time periods 
T a ST permissible delivery window for unit u with aircraft a and route 
r combination; this time window covers the period from the unit 


is Available-to-Load Date to the Required Delivery Date plus a 
maximum allowed lateness for delivery 


3. Data 


a. Movement Requirement Data 


MovePAX Troop movement requirement for unit u from origin i to 
destination K 
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MoveUE,,, 


ProBulk, 


ProOver, 


ProOut, 


Equipment movement requirement in tons for unit 4 from origin 
i to destination K 


Proportion of unit 4 cargo that is bulk-sized (loose cargo 
palletized on a 88" x 108" platform) 


Proportion of unit u cargo that is over-sized (non-palletized cargo 
rolling stock that can fit into a C141); over-sized cargo has a 
larger dimension than bulk-sized cargo 


Proportion of unit u cargo that is out-sized (non-palletized cargo 
that can fit into a C5 or C17 but not a C141) 


b. | Penalty Data 


LatePenUE, 
LatePenPAX, 
NoGoPenUE, 
NoGoPenPAX, 


MaxLate 


Lateness penalty (per ton per day) for unit u equipment 
Lateness penalty (per soldier per day) for unit u troops 
"Non-delivery" penalty (per ton) for unit u equipment 
"Non-delivery" penalty (per soldier) for unit u troops 


Maximum allowed lateness (in days) for delivery 


c. Cargo Data 


UESqFt, 


PAXWt 


Average cargo floor space (in sq. ft.) per ton of unit u equipment 


Average weight of a soldier inclusive of personal equipment 


d. Aircraft Data 


Supply, 
MaxPAX, 


PAXSgFt, 


ACSqFt, 


Additional number of aircraft a made available on day t 
Maximum troop carriage capacity of aircraft a 


Average cargo space (in sq. ft.) consumed by a soldier for 
aircraft a 


Cargo floor space (in sq. ft.) of aircraft a 
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LoadEff, Cargo space loading efficiency for aircraft a 


URate, Established utilization rate (flying hours per aircraft per day) for 
aircraft a 
EffTime, Number of effective time periods for an aircraft that is available 


from time period 1; EffTime,=T-t+1 where T is the time horizon 


e. Airfield Data 


MOGCap, "Maximum On Ground" capacity for narrow-body aircraft at 
airfield b 
MOGReq,, Normalized MOG requirement (to narrow-body aircraft) for 


aircraft a at airfield b. Normalization is required for comparison 
purposes since an airfield has different MOG figures for each 
aircraft type (body type); i.e., an airfield can accommodate 
different number of each type of aircraft (or combinations) 


MOGEff MOG efficiency factor; introduced as a discount factor as it is 


impossible to fully utilize the MOG capacity due to the stochastic 
nature of aircraft ground time 


f. Aircraft Route Performance Data 
MaxLoad,, Maximum payload (in tons) for aircraft a flying route r. 


DGTime,,, Aircraft ground time (due to onload of cargo, enroute refuelling, 
maintenance, etc) for aircraft a at airfield b on delivery route r 


RGTime,,, Aircraft ground time taken by aircraft a at airfield b on recovery 
route r 
DTime,,, Cumulative time (flight time plus ground time) taken by aircraft 


a to reach airfield b along delivery route r 


RTime,,, Cumulative time (flight time plus ground time) taken by aircraft 
a to reach airfield b along recovery route r 


FltTime,, Total flying hours consumed by aircraft a on route r 
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DaysLate,,, Number of days late (after required delivery date) in delivering 
unit 4 requirement using aircraft a flying route r with mission 
start time z; DaysLate,,, is equal to zero if the cargo/troops 
arrive at the destination before (or on) the required delivery date 


4. Decision Variables 


Xuan Number of aircraft a that airlift unit u, via route r with mission 
start time / 
Y... Number of aircraft a that recover from a destination airfield via 


route r with start time { 


Allot,;, Additional number of aircraft a available from day t that are 
allocated to origin / 


Has Number of aircraft a inventoried at origin i, at time / 
IB. Number of aircraft a inventoried at destination K, at time / 
TonsUE,, Total tonnage of unit 4 equipment airlifted by aircraft a, via 


route r with mission start time t 


TPAX an Total number of unit u troops airlifted by aircraft a, via route r 
with mission start time / 


UENoGo,,, Total tonnage of unit 4 equipment not airlifted from origin i to 
destination K in the prescribed time frame 


PAXNoGo,; Number of unit u troops that not airlifted from origin 7 to 
destination K in the prescribed time frame 


5. Objective Function 
The model objective is to minimize the total weighted penalties (weighted by 
the movement priority of a unit) for late deliveries and undelivered cargo and troops. 


The objective function value, z, is given by: 
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Minimize Z = 


Y Y Y Y (LatePenUE, * DaysLate,,,,) * TonsUE,;, 


u a reR, IET 


+ Y Y Y Y, (atePenPAX, * DaysLate,,,) * TPAX,,, 


u a FER TET 


n y ». > (NoGoPenUE_ * UENoGo,,,) + (NoGoPenPAX, * PAXNoGo,,). 


u I 


There are two reasons for including a "undelivered" cargo category in the 
objective function. First, it allow the user to control the number of time periods. 
Second, it permits situations where the movement requirements cannot be fulfilled due 
to limited aircraft, airfields, and time horizon; i.e., it allows the LP model to produce 
a feasible solution even when available assets are inadequate. 

Some care must be taken in determining late and non-delivery penalties to 
ensure that they are on a common basis of x tons of one unit versus x tons of another. 
In addition, the user should ensure that weights are consistent; i.e., late delivery should 
be preferred to non-delivery. The weights will be consistent provided a late penalty (per 
ton per day) is less the corresponding non-delivery penalty (per ton) divided by the 


maximum allowed lateness (in days) for delivery. 
6. Constraints 


The objective function is subject to the following constraints. 


a. Demand Satisfaction Constraints 
Demands on the airlift system are the movement requirements for troop 


and different classes of cargo (bulk, over and out-sized). There are four different kinds 
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of demand satisfaction constraints (one for troops and three for cargo). Three separate 
constraints are required for cargo due to the need to impose upper bounds on the amount 
of cargo that can be carried by a particular class of carrier (e.g., over-sized cargo 
carrier); this ensures cargo-carrier compatibility. 

The demand satisfaction constraints for troops ensure that the number 
of troops airlifted within the permitted time window for each unit is greater than or equal 
to the number required, with the shortfalls accounted for in the deviation variables 
(PAXNoGo). The demand satisfaction constraints for each unit's cargo serve the same 


purpose for the required equipment. The demand satisfaction constraints follow: 


(1) Demand Satisfaction Constraints for All Classes of Cargo. For 
each unit and associated origin-destination pair, the total tonnage of unit equipment 
delivered to the theater by all aircraft types (over all permissible aircraft, route and 
mission start time combination) plus total tonnage not delivered (UENoGo) must be 


greater than or equal to the movement requirement for equipment: 


y Y y TonsUE,,,, * UENoGo,, 2 Μονεῦξιι ο CICER: MoveUE, , ? 0. 


αι. τες A 

(2) Demand Satisfaction Constraints for Out-Sized Cargo. For each 
unit and associated origin-destination pair, the total tonnage of unit equipment carried by 
out-sized cargo capable aircraft plus the total tonnage not delivered must be greater than 
or equal to the movement requirement for out-sized cargo (given by proportion of out- 


sized cargo multiplied by movement requirement for equipment): 
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y y y TonsUE,,, * UENoGo,, 2 ProOut, * MoveUE,,, 


a cA rem te 1227 
V u,ik : MoveUE, , ? O. 
(3) Demand Satisfaction Constraints for Over-Sized Cargo. For each 
unit and associated origin-destination pair, the total tonnage of unit equipment carried by 
over-sized cargo capable aircraft plus the total tonnage not delivered must be greater than 


or equal to the movement requirement for over-sized and out-sized cargo: 


Y Y Y, TonsUE,,, * UENoGo,, > (ProOver, * ProOut,) * MoveUE,,, 


acA,, TER 16 1, 
V u,i,k : MoveUE, , >0. 
The right-hand-side of the constraint is valid because out-sized 


cargo capable carriers are also over-sized cargo capable. 


(4) Demand Satisfaction Constraints for Troops. For each unit and 
associated origin-destination pair, the total number of troops airlifted to the theater plus 


those not airlifted must be greater than or equal to the movement requirement for troops: 


Σ y Y TPAX,, + PAXNoGo,, = MovePAX,, — V uj: MovePAX,, 0. 


α ΤΕΝ, τε]. 
b. Aircraft Balance Constraints 
The number of aircraft a assigned for airlift missions (or to be 
recovered from a destination airfield) must be no more than the number of aircraft a 


avallable for use. The three different kinds of aircraft balance constraints are: 
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(1) Aircraft Balance Constraint at Origin Airfield. For each aircraft 
a, origin i and time period ¢ combination, the total number of aircraft assigned for airlift 
missions, plus those inventoried for later use must be equal to the total number of aircraft 
available from the previous period plus new supply of aircraft allocated to the origin and 


returns from previous missions: 


23 > Αρ y H E H ia-1) p Allot 5 y > De V. aj. 


u r €DR, reR, t! + RTime ,,, = t 


Note: H,,,_) is only defined for t greater than 1. 


X ar, 35 Only defined for teT,,,. 


uart 

(2) Balance Constraint for Aircraft Allocation. The total number of 
aircraft a allocated to the different origin airfields in each time period ¢ must not be 
greater than the new supply of aircraft a. This constraint is incorporated in the model 
as it is more efficient for the LP model to distribute the new supply of aircraft than for 


a user to predetermine the allocation: 


y Allot; < Supply, V ay : Supply, » 0. 


(3) Aircraft Balance Constraint at Destination Airfield. For each 
aircraft a, destination k and time period ¢ combination, the total number of aircraft a 
recovering from the destination airfield this period plus those that are to be inventoried 
must be equal to those waiting to return (inventory from last period) plus new arrivals 


at the destination: 
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Y Yn + HPA == HP  »ὁ ο Ἱ -ι 


re RR, u r€R, t'ET y 


t'+DTime „=t 


Note: HP akte-1) 


is only defined for t greater than 1. 
c. X Aircraft Physical Limitation Constraints 
There are three different kinds of constraints on the physical limitations 


of aircraft. These constraints ensure that limits on troop carriage capacity, maximum 


payload, and cargo floor space are observed. 


(1) Troop Carriage Capacity Constraints. The total number of troops 
airlifted for each unit u, aircraft a, route r, and mission start time t combination must not 
be greater than the troop carriage capacity of aircraft a multiplied by the number of 


aircraft a assigned: 


TPAX , Ss MaxPAX, *X_, V u,a,r,t : t€ T 


rt uar° 


(2) Maximum Payload Constraints. The total payload (equipment and 
troops) airlifted for each unit u, aircraft a, route r, and mission start time t combination 
must not be greater than the maximum payload for aircraft a flying route r, multiplied 


by the number of aircraft a assigned: 


TonsUE,,, * (PAXWt * TPAX, ,) «s MaxLoad, * X, V uart: c 


rt 


(3) Cargo Floor Space Constraints. The total floor space taken up 


by troop and equipment for each unit u, aircraft a, route r, and mission start time { 
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combination must not be greater than the aircraft cargo floor space multiplied by the 
loading efficiency factor and number of aircraft a assigned: 
(PAXSqFt, * TPAX, ,) + (UESqFt, * TonsUE,,,) < ACSqFt, * LoadEff, * X H 


V UOI ved os 


d. Aircraft Utilization Rate Constraints 

These constraints ensure that the average flying hours consumed per 
aircraft per day is less than AMC's established utilization rate for each aircraft type; the 
upper bounds for aircraft utilization is designed to capture spares availability, aircraft 
reliability, crew availability, etc. 

Aircraft utilization rate constraints are modelled by comparing the actual 
flying hours consumed by an aircraft fleet in a deployment with the maximum achievable 
flying hours for the fleet based on the utilization rate. The total flying hours consumed 
by an aircraft fleet in a deployment is equal to the total flying hours consumed for both 
delivery and recovery routes. The maximum achievable flying hours is given by the 


utilization rate multiplied by the productive time periods of each available aircraft. 


Y Y Y (FiTim, * X4) + 9», 5  (FltTime, + Y,,,) 


u ren; ref reR, t 


< Ë. URate, * Efflime, * Supply,, V a. 
t: Supply ;,> 0 


As a simple illustration for the above equation, consider a fleet of 5 
aircraft made available from day 11. If the utilization rate for this aircraft type is 10 


flying hours per aircraft per day and the time window of concern is 30 days, then the 
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maximum achievable flying hours for the aircraft type over the time period is 1000 (10 
x 20x 5). This total must not be exceeded by the actual flying hours consumed by the 


aircraft fleet over the same period. 


e. Aircraft Handling Capacity of Airfield (MOG Constraint) 

These constraints model the aircraft handling capacity or throughput of 
each airfield. They ensure that the LP does not route more aircraft through an airfield 
than it can handle each day. 

The maximum number of aircraft an airfield can handle each day 
depends on the following factors: MOG capacity of the airfield which gives the 
maximum number of aircraft it can accommodate at any one time, aircraft body type 
(narrow-body, wide-body, etc) and the amount of time each aircraft spends at the 
airfield. 

Aircraft handling capacity constraints are modelled by computing the 
total MOG consumed (normalized to that of a narrow-body aircraft) in each period by 
each individual aircraft using the airfield. This total should not exceed the normalized 
MOG capacity of the airfield discounted by a MOG efficiency factor. The MOG 
efficiency factor is introduced as the stochastic nature of aircraft ground time implies that 
it is impossible to fully utilize the airfield. It is recognized that deviations from 
scheduled ground times can cause some aircraft to wait on the ground (or conceivably 
in the air) while others are being serviced. The MOG consumed by an aircraft is equal 
to its normalized MOG requirement multiplied by the ratio of time spent (in hours) under 


service to number of hours in a time period. The equation follows: 
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y y y 3. (MOGReq,, * DGTime,, | 24) * X, 
u a Lem. "OT. 


t «DTime,,, =f 


ΕΤ τ dE (MOGReg,, * RGTime,, | 24) * Y. 


a reR, 1'»RTime,, =1 
< MOGEff * MOGCap, V bt. 

As a simple illustration for the MOG constraint, consider the case of 
an airfield that can handle 10 narrow-body aircraft at a time and a MOG efficiency factor 
of 0.8. If each narrow-body aircraft spends 3 hours on the ground, then over the whole 
period (day), the airfield would be able to handle 64 narrow-body aircraft; 1.e., as the 
equation implies 1 x 3/24 x 64 — 0.8 x 10. 

An assumption made in the above formulation is that aircraft 
inventoried at origin or destination airfields do not consume any MOG capacity. This 
IS not strictly valid since inventoried aircraft take up parking space. However, since 
MOG figures are based on a host of other factors beside parking spaces, it is probably 
more accurate to model the problem as such. As noted in the Assumptions section in 
Chapter II, if the USAF could provide separate data for "parking space MOG" and 


"ground services MOG", a more accurate modelling approach would be possible. 


C. SUMMARY 
The mathematical expressions for the objective function and constraints are 


summarized below: 
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€ Objective Function 


Minimize Z = 
Y > Y Y (LatePenvE, * DaysLate,,.) * TonsUE,,,, 


u a r€ R, tr 


+ Y Y Y Y GatePenPAX, * DaysLate,,,) * TPAX,,, 


a reR, led 


$3099 » > (NoGoPenUE, * UENoGo,,,) + (NoGoPenPAX, * PAXNoGo, ,). 


€ Constraints 


Y” Y Y TonsUE,, * UENoGo,, 2 MoveUE,, V  u,i,k : MoveUE,,,>0 


σΕΑ, τον τει, 


Y” Y Y, TonsUE,,, * UENoGo,, = ProOut, * MoveUE,, 
QEA u TER tel. 


V u,t,k : MoveUE_ >0 


Y); Y Y TonsUE,,, + UENoGo,, » (ProOver, * ProOut,) * MoveUE,,, 
aCA TERI TET 


y u,ik : Move UE 20 


Y ` Y TPAX, * PAXNoGo,, > MovePAX,,  W u,i,k : MovePAX,,>0 


a reKR,, t€ T.z, 


Y Y X. * H; ° Hanlon, κ» PED V Xl 


uart ait 
u reDR, rER, 1 +RTime,, =1 
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Y” Allot, « Supply, V a: Supply, > 0 


πο ο πο UMP E c EE Dy 6 gas 


r€ RR, u reR, Τετ, 
t! « DTime,, 7t 
TPAX 4, $ MaxPAX * X,, V Halt cp 
TonsUE, , * (PAXWt * TPAX, ,) < MaxLoad, * X, V ματι tel 


(PAXSqFt, * TPAX,,,) + (UESqFt, * TonsUE,,,) s ACSqFt, * LoadEff, * X, 


V αι tel. 


ο ο M μπε Χο + 5 V (Fifime, + Y) 


u reR, eT FER, t 


< P URate, * Efflime, * Supply, V a 
t: Supply „> 0 


y 9 y M (MOGReq,, * DGTime,,, | 24) * X, 
u a reR, εἰ, 


t'+DTime,,,=t 


+ y 35 y (MOGReg,, * RGTime,,, | 24) + Y, 


a TER, i('«RTime,, 7t 


« MOGEff * MOGCap, V bt 


X ,Y .Allot 


uart? an? ait? 


H 


ait? 


HP, TonsUE, ,, TPAX,, ,, UENoGo 


uart? uart? uik? 


PAXNoGo,, > 0 


V arp 
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IV. IMPLEMENTATION DETAILS 
This chapter discusses the General Algebraic Modelling System (GAMS) 
implementation of the conceptual optimization model described in Chapter III and can 
be skipped without any loss in continuity. The materials of this chapter include key 


implementation features, and a mathematical formulation of the implemented model. 


A. ROUTE IMPLEMENTATION AND INDEX SUBSETS 
The data structures used to implement air route structure and index subsets facilitate 


minimal external data processing. 


]. Route Implementation 
The air route structure of the optimization model is represented by five 
indices as opposed to a single index as described in the conceptual model. These five 
indices represent the origin, destination and three enroute airfields. The use of five 
airfield indices enable the model to distinguish the origin, enroute, and destination 
airfields and thus allows aircraft flight times and ground times to be calculated more 
easily. A graphical illustration of the air route structure is shown in Figure 1 (characters 


in the circular nodes are airfield GEOLOC codes). 
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Origin Enroute 1 Enroute 2 Enroute 3 Destination 
Airfields Airfields Airfields Airfields Airfields 
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Figure 1. Air Route Structure 


Although three indices are used to denote enroute stops, this does not imply 
that a delivery or recovery route must have three enroute stops. In fact, the model 
supports delivery/recovery routes with zero to three enroute stops. Routes with less than 
three enroute stops are implemented by replacing the non-existent enroute airfield names 
with the destination airfield name; e.g., the delivery route shown in Figure 1 has two 


enroute stops. 


2. Index Subsets 
Some of the index subsets of the conceptual model described in the previous 
chapter are not explicitly declared as a subset in the implemented model; i.e., they are 
not index subsets per se. Instead, these virtual index subsets are derived from the main 


index set through logical conditional checks performed by GAMS in model generation. 
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For example, the permitted delivery time window for each unit (a virtual subset) is 
generated by the model using relevant data such as the unit available load date, the 
required delivery date, and the time taken by an aircraft type to deliver the cargo. The 
generation of these induced virtual index subsets will become apparent in the objective 


function and constraint equations described below. 


B. MATHEMATICAL FORMULATION 

The mathematical formulation of the optimization model as implemented in GAMS 
is given below. To facilitate a concise discussion, only index and data sets explicit in 
the objective function or constraints will be described. The interested reader can refer 


to the GAMS formulation in Appendix B for additional detail. 


1. Indices 


u unit; e.g., 82nd airborne 

a alrcraft type; e.g., C5, C141 

ab aircraft body type; e.g., wide-body (wb), narrow-body (nb) 

ga ground activity; e.g., load (onld), off-load (offld), enroute (enr) 
turnaround 

t,tp time period in days 

c cargo type; e.g., unit equipment (UE), troops (PAX) 

cc cargo class; e.g., bulk, over or out-sized cargo 

af,afp generic airfield 

i origin airfield; IS AF 

k destination airfield in the theater; KS AF 

el first enroute airfield; EI ΑΕ 

τ second enroute airfield; E2C AF 

e3 third enroute airfield; E3 © AF 


Note: A compound index r will be used as a compact notation to represent the 
airfield indices (1,k,e1,e2,e3) in later discussions whenever appropriate. 
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2. Data 


a. Movement Requirement Data 


MovePAX(u,i,k) 


MoveUE(u,i,k) 


CargoP(u,cc) 
ALD(u) 


RDD(u) 


Troop movement requirement (in 100s) for unit u from 
origin í to destination K 


Equipment movement requirement (in 100 tons) for unit u 
from origin í to destination k 


Proportion of cargo classes (bulk, over, out) of unit u 
Available load date for unit u 


Required delivery date for unit u 


b. Penalty Data 


LatePen(u,c) 


NoGoPen(u,c) 


MaxLate 


c. Cargo Data 


UESqFt(u) 


PAX Wt 


Lateness penalty for cargo type c of unit u. LatePen for 
troops is measured in terms of per 100 PAX per day and 
for equipment is in terms of per 100 ton per day. 


Non-delivery penalty for cargo type c of unit u. NoGoPen 
for troops is measured in terms of per 100 PAX and for 
equipment is in terms of per 100 ton. 


Maximum allowed lateness for delivery. This maximum 
lateness is necessary both for controlling the size of the 
time domain and to allow the LP to produce a feasible 
solution even when assets are inadequate. 


Average cargo space (in 1000 sq. ft.) taken up by 100 tons 
of unit u’s equipment 


Average weight (in 100 tons) of 100 troops inclusive of 
personal equipment 
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d. Aircraft Data 


Supply(a,t) 


ACCargo(a,cc) 


MaxPAX(a) 


PAXSgFt(a) 


ACSgqFt(a) 
LoadEff(a) 


GTime(a,ga) 


URate(a) 


Additional number of aircraft a made available on day t 


Aircraft-cargo class matching table; a value of 1) 
indicates compatibility 


Maximum troop carriage capacity (in 100s) of aircraft a 


Average cargo space (in 1000 sq. ft.) consumed by 100 
soldiers for aircraft a 


Cargo floor space (in 1000 sq. ft.) of aircraft a 
Cargo space loading efficiency for aircraft a 


Ground time required to accomplish ground activity ga for 
aircraft a 


Established utilization rate (in 100 flying hours per aircraft 
per day) for aircraft a 


e. Airfield Data 


MOGCap(af,"nb") 


MOGReq(af,a) 


MOGEff 


Dist(af,afp) 


"Maximum On Ground" capacity for narrow-body aircraft 
at airfield af 


Normalized MOG requirement (to narrow-body aircraft) 
for aircraft a at airfield af 


MOG efficiency factor; introduced as a discount factor as 
it is impossible to fully utilize the MOG capacity due to 


the stochastic nature of aircraft ground time 


Distance (in nautical miles) between airfields af and afp 


f. Aircraft Route Performance Data 


MaxLoad(a,r) 


Maximum payload (in 100 tons) for aircraft a flying on 
route r 
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CTtoE1(a,r) 


CTtoE2(a,r) 


CTtoE3(a,r) 


CTtoK(a,r) 


RCTtoE3(a,r) 


RCTtoE2(a,r) 


RCTtoEl(a,r) 


RCTtol(a,r) 


FltTime(a,r) 


VRouteX(a,r) 


VRouteY (a,r) 


Cumulative time (in hrs) taken by aircraft a flying on 
delivery route r to reach enroute e7. This cumulative time 
includes both flying time and ground time. 


Cumulative time (in hrs) taken by aircraft a flying on 
delivery route r to reach enroute e2 


Cumulative time (in hrs) taken by aircraft a flying on 
delivery route r to reach enroute e3 


Cumulative time (in hrs) taken by aircraft a flying on 
delivery route r to reach destination airfield k 


Cumulative time (in hrs) taken by aircraft a flying on 
recovery route 7 to reach enroute e3 


Cumulative time (in hrs) taken by aircraft a flying on 
recovery route r to reach enroute e2 


Cumulative time (in hrs) taken by aircraft a flying on 
recovery route r to reach enroute el 


Cumulative time (in hrs) taken by aircraft a flying on 
recovery route r to reach origin i 


Total flight time (in 100 hrs) consumed by aircraft a flying 
on route r; excludes ground time 


A dynamic set with members indicating aircraft-delivery 
route compatibility. This set is used to control aircraft 
routing; e.g., military aircraft fly military routes and 
civilian aircraft fly civilian routes. 


A dynamic set with members indicating aircraft-recovery 
route compatibility 


3. | Decision Variables 


X(u,a,r,t) 


Number of aircraft a assigned to airlift unit u, via route r 
with mission start time z. An aircraft is permitted to carry 
only a single unit's cargo and troops as in a typical airlift 
mission 
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Y(a,r,t) 


Allot(a,1,t) 


H(a,i,t) 
HP(a,k,t) 


TonsUE(u,a,r,t) 


TPAX(u,a,r,t) 


UENoGo(u,1,k) 


PAXNoGo(u,i,k) 


4. Dynamic Sets 


Number of aircraft a that recover from a destination 
airfield via route r with start time { 


Additional number of aircraft a available from day t that 
are allocated to origin i 


Number of aircraft a inventoried at origin i at time { 
Number of aircraft a inventoried at destination k at time £ 


Total tonnage (in 100 tons) of unit u equipment airlifted by 
aircraft a, via route r with mission start time t 


Total number of unit 4 troops airlifted by aircraft a via 
route r with mission start time t 


Total tonnage (in 100 tons) of unit 4 equipment not 
airlifted from origin i to destination K in the prescribed 
time frame 


Number of unit u troops (in 100s) not airlifted from origin 
i to destination K in the prescribed time frame 


Two dynamic sets, DSetX(u,a,r,t) and DSetY(a,r,t), are used repeatedly 


throughout the formulation of the objective and constraint equations to reduce the number 


of decision variables and thus the size of the LP. 


The dynamic sets serve two purposes. First, they control the allowable 


combination of unit u, aircraft a, route r and mission start time ¢ for airlift missions. 


Second, the dynamic sets eradicate certain combinations of decision variables that are 


known to have zero values at optimality; 1.e., instead of keeping these decision variables 


in the LP model and letting the optimizer drive them to zero, these decision variables are 


removed from the model formulation. 
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a. DSetX(u,a,r,t) 


DSetX(u,a,r,t) is used to control the valid combination of unit u, 


aircraft a, route r and mission start time ¢ for the following decision variables: 


X(u,a,r,t), TonsUE(u,a,r,t) and TPAX(u,a,r,t). For a combination to be valid, it must 


satisfy the following conditions: 


Mission start time t must satisfy the following two inequalities: t 2 ALD(u) and 
t € RDD(u) - [CTToK(a,r)/24] 4 MaxLate; i.e., an unit can only be airlifted if 
it is ready to move and that the cargo will reach its destination before the required 
delivery date plus the maximum allowed lateness. Square brackets denote a "round 
operator"; e.g., [3.7] = 4 and [3.4] = 3 


There is movement requirement for unit u from origin airfield í to destination 
airfield k; i.e., MoveUE(u,i,k) > O or MovePAX(u,i,k) > 0 


Be a member of VRouteX(a,r); i.e., permitted aircraft-delivery route combination 
A particular aircraft-route combination is to be considered only if the permissible 
payload for the aircraft is at least 25 percent of the aircraft full load capability. 
This 25 percent factor is for cost efficiency and is also representative of real airlift 


operations. 


Supply of aircraft type up to time t 15 greater than zero 


b. DSetY(a,r,t) 


DSetY (a,r,t) is used to control the valid combination of aircraft a, route 


r and start time í for the recovery decision variable Y(a,r,t). For a combination to be 


valid, it must satisfy the following conditions: 


e Be a member of VRoute Y (a,r); 1.e., permitted aircraft-recovery route combination 


e Supply of aircraft type up to time 1 is greater than zero 


33 


5. Objective Function 
The model objective is to minimize the total weighted penalties (weighted by 
the movement priority of a unit) for late deliveries and undelivered cargo and troops. 
The objective function, z, is given by: 
Minimize Z = 


ΣΣΣ 


r t€DSetX(u,a,rX) 


( LatePen(u," UE") x Max(0, (t * [CTtoK(a,r|24] - RDD(u))) * TonsUE(u,a,r,t) 
* (LatePen(u," PAX ^ x Max(0, (t - [CTtoK(a,r)|24] - RDD(u))) * TPAX(u,a,r,t) 


ΣΣ 


u i k:MoveUE(u,;k)>0orMovePAX(u,i,k)>0 


(NoGoPen(u, UE") » UENoGo(u,i,k)) + (NoGoPen(u,”"PAX”) x« PAXNoGo(u,i,k)). 


Note. 2:200, Degli 


k el e2 e3 


The number of days late in delivery for a particular unit u, aircraft a, route 
r and mission start time t combination Is given by the GAMS function Max(a,b) where 
a and b are numerical values or arithmetic expressions. For cargo and troops delivered 
to the theater on or before the required delivery date no penalty is incurred. For cargo 
and troops delivered late, the number of days late in delivery Is equal to the mission start 
time plus time taken by the aircraft to deliver the cargo and troops minus the required 


delivery date. 


6. Constraints 


The objective function is subject to the following constraints. 
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a. Demand Satisfaction Constraints 
The constraints for satisfying demand are the requirements for troops 
and the different classes of cargo (bulk, over and out-sized). There are four different 
kinds of demand satisfaction constraints (one for troops and three for cargo). Three 
separate constraints are required for cargo due to the need to ensure aircraft-cargo 


compatibility. The demand satisfaction constraints follow: 


(1) Demand Satisfaction Constraints for All Classes of Cargo. For 
each unit and associated origin-destination pair, the total tonnage of unit equipment 
delivered to the theater by all aircraft types (over all permissible route and mission start 
time combination) plus total tonnage not delivered (UENoGo) must be greater than or 


equal to the movement requirement for equipment: 


y O Y TonsUE(u,a,r,t) + UENoGo(u,i,k) 


a:ACCargo(a,“bulk”)>0 εἰ ε εὖ teDSetX(ua,r,t) 
> MoveUE(u,i,k) V uk : MoveUE(u,i, k) » 0. 

(2 Demand Satisfaction Constraints for Out-Sized Cargo. For each 
unit and associated origin-destination pair, the total tonnage of unit equipment carried by 
out-sized cargo capable aircraft plus the total tonnage not delivered must be greater than 
or equal to the movement requirement for out-sized cargo (given by proportion of out- 


sized cargo multiplied by movement requirement for equipment): 


p» y y y y TonsUE(u,a,r,t) + UENoGo(u,i,k) 


a:ACCargo(a,"out”)>0 el el eJ t€ DSetX(u,a,r,t) 


> CargoP(u,“out”) * MoveUE(u,i,k) V u,ik : MoveUE(u,i,k) > 0. 
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(3) Demand Satisfaction Constraints for Over-Sized Cargo. For each 
unit and associated origin-destination pair, the total tonnage of unit equipment carried by 
over-sized cargo capable aircraft plus the total tonnage not delivered must be greater than 


or equal to the movement requirement for over and out-sized cargo: 


» ΣΣΣ Y  TonUEG,,) * UENoGo(ui,O 


a:ACCargo(a, "over ) » 0 el e2 e3 t E DSetX(u,a,r,t) 
> (CargoP(u,”over") + CargoP(u,"out")) * MoveUE(u,i,k) 
V u,,k : MoveUE(u,i,k) >0. 
The right-hand-side of the constraint is valid because out-sized 


cargo capable carriers are also over-sized cargo capable. 


(4) Demand Satisfaction Constraints for Troops. For each unit and 
associated origin-destination pair, the total number of troops airlifted to the theater plus 


those not airlifted must be greater than or equal to the movement requirement for troops: 


3 y y y a TPAX(u,a,r,t) + PAXNoGolu,i,k) > MovePAX(u,t,k) 


a el e2 e3 teDSetX(u,a,r,t) 
V u,,k : MovePAX(u,i,k)>0. 
b. Aircraft Balance Constraints 
The number of aircraft a assigned for airlift missions (or to be 
recovered) must be no more than the number of aircraft a available for use. The three 


different kinds of aircraft balance constraints are: 


(1) Aircraft Balance Constraint at Origin Airfield. For each aircraft 


a, origin i and time period ¢ combination, the total number of aircraft assigned for airlift 
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missions, plus those inventoried must be equal to the total number of aircraft available 
from the previous period plus new supply of aircraft allocated to the origin and returns 


from previous missions: 


y » y y y ματι + H(a,i,t) 
u k 


el e2 e3€ DSetX(u,a,r,t) 


= H(a,i,t-1) + Allot(a,i,t) + y Y Y >` y Y(a,r ip) VW αμ. 
k 


el e? e3 tp € DSetY(a,r,tp) 
tp + [RCTtol(a,r)/24] =t 


Note: H(a,i,t-1) is only defined for t greater than 1. 
(2) Balance Constraint for Atrcraft Allocation. The total number of 
aircraft a allocated to the different origin airfields in each time period t must be no more 


than the new supply of aircraft a: 


Y” Allot(ai,)) < Supply(a,t) V — ay : Supply(a,t) »O. 


(3) Aircraft Balance Constraint at Destination Airfleld. For each 
aircraft a, destination k and time period t combination, the total number of aircraft a 
recovering from the destination airfield this period plus those inventoried must be equal 
to those waiting to return (inventory from last period) plus new arrivals at the 


destination: 


»» > > y Y(a,r,t) + HP(a,k,t) 


l el e2 e3eDSetY(a,r,t) 


- HP(akt-1 + Y Y Y Y y` X(u, ar ip) Y akt. 


l el e2 ελ tp € DSetX(u,a,r,tp) 
tp * [CTtok(a,r)24] * t 


Note: HP(a,k,t-1) is only defined for t greater than 1. 
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c. Aircraft Physical Limitation Constraints 
The aircraft physical limitation constraints modelled are the troop 


carriage capacity, maximum payload and cargo floor space constraints. 


(1) Troop Carriage Capacity Constraints. The total number of troops 
airlifted for each unit u, aircraft a, route r, and mission start time £ combination must not 
be greater than the troop carriage capacity of aircraft a multiplied by the number of 
aircraft a assigned: 

TPAX(u,a,r,t) < MaxPAX(a) * X(u,a,r,t) V  u,a,r te DSetX(u,a,r,t). 

(2) Maximum Payload Constraints. The total payload (equipment and 
men) airlifted for each unit u, aircraft a, route r, and mission start time £ combination 
must not be greater than the maximum payload for aircraft a flying route r, multiplied 
by the number of aircraft a assigned: 

TonsUE(u,a,r,t) + (PAXWt * TPAX(u,a,r,t)) < MaxLoad(a,r) * X(u,a,r,t) 
V u,ya,r,t € DSetX(u,a,r,t). 

(3) Cargo Floor Space Constraints. The total floor space taken up 
by troops and equipment for each unit uw, aircraft a, route r, and mission start time t 
combination must not be greater than the aircraft cargo floor space multiplied by the 
loading efficiency factor and number of aircraft a assigned: 

(PAXSqFt(a) * TPAX(u,a,r,t)) + (UESqFt(u) * TonsUE(u,a,r,t)) 


< ACSqFt(a) * LoadEfRa) * X(u,a,r,t) 
V uar, €DSetX(u,a,r,)), ACSqFt(a)»0. 
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d. | Aircraft Utilization Rate Constraints 

These constraints ensure that the average flying hours consumed per 
aircraft per day is less than AMC’s established utilization rate for each aircraft type. 
Aircraft utilization rate constraints are modelled by comparing the actual flying hours 
consumed by an aircraft fleet in a deployment with the maximum achievable flying hours 
based on the utilization rate. The total flying hours consumed by an aircraft fleet in a 
deployment is equal to the total flying hours consumed on both delivery and recovery 
routes. The maximum achievable flying hours for each aircraft is given by the utilization 
rate multiplied by the number of productive time periods of that aircraft (CARD(t) 


represent the time horizon in the equation). 


y y y FliTimela,r) * X(u,a,r,t) + y y FltTimela,r) * Y(a,r,t) 


u r t€DSetX(u,a,r,t) r 1€DSetY(a,r,t) 


< y URate(a) * CCARD(t) - 1 -t) *Supply(at) V a. 


t: Supply(a,t) » O 


e. Aircraft Handling Capacity of Airfield (MOG Constraint) 

These constraints model the throughput or handling capacity of each 
airfield. They ensure that the LP does not route more aircraft through an airfield than 
it can handle each day. 

The aircraft handling capacity constraints are modelled by computing 
the total MOG consumed (normalized to that of a narrow-body aircraft for comparison) 
in each period by each individual aircraft using the airfield. An airfield can serve as an 
origin, enroute, or a destination airfield; thus the total MOG consumed at a particular 


airfield on each day is the sum of the individual consumptions when serving in these 
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different capacities. The total MOG consumed must be less than the normalized MOG 
capacity of the airfield discounted by the MOG efficiency factor (it is impossible to fully 
utilize the MOG capacity due to deviations from scheduled ground times). 

The airfield capacity that an aircraft consumes is equal to its MOG 
requirement multiplied by the ratio of time spent (in hours) under service to number of 
hours in a time period (in this case, 24 hours). MOG requirement are normalized to that 
of a narrow-body aircraft. The amount of time an aircraft spends at origin, enroute and 
destination airfields are the loading time (onload), enroute refuelling time (enr) and off- 
load time (offld) respectively. In the formulation of the constraints, MOG consumption 
at the origin and destination airfields are associated with the X(u,a,r,t) variables only 
whereas, the MOG consumption at enroute airfields are associated with both the 
X(u,a,r,t) and Y(u,a,r,t) variables. 

As the model supports air routes with zero to three enroute stops, it is 
also necessary to ensure that there is no double counting of MOG consumption 
(destination airfield names are used to replace non-existent enroute airfield names). 
Potential double counting of MOG consumption is prevented through conditional checks 


performed by the constraint equations. The constraint equation follows: 
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+ 


+ 


+ 


+ 


22222 


k εἰ e2 e3eDSeiX(u,a,af k,el,e2,e3,t) 
(MOGRegq(af,a) « GTime(a,”onld")/24) * X(u,a,af,k,e1,e2,e3,t) 


ΣΣ 2 


e2 e3 ipeDSeiX(u,a,ik,af,e2,e3,tp) 
tp * [CTtoEl(a,r)/24) *t 
Dist(i,af) 0 
CTtoEl(a,r) » CTtoK(a,) 


(MOGReq(af,a) * GTime(a,"enr") /24) * X(u,a,i,k,af,e2,e3,tp) 


2222) 


e2 e3  tpeDSetY(a,ik af,e2,e3,tp) 
tp * [RCTtoE1(a,r)/24] *t 
Dist(af,e2) + O 
RCTYoEl(a,y) » RCTYol(a,r) 


( MOGReq(af,a) * GTime(a, "enr "^)/24) * Y(a,i,k,af,e2,e3,tp) 


222222 


k el e3 tpeDSetX(u,a,i.k.el af,e3 tp) 
tp + [CTtoE2(a,r)/24] =1 
Dist(el,af) * 0 
CTtoE2(a,r) « CTtoK(a,r) 


( MOGReq(af,a) * GTime(a," enr ")/24) » X(ua,i,k,e1,af,e3,tp) 


222.22 > 


k el e3 tpeDsSetY(aak,el,af,e3,tp) 
tp * [RCTtoE2(a,r)/24] * t 
Dist(af,ee3) * O 
RCTioE2(a,r) += RCTiol(a,r) 


( MOGReq(af,a) » GTime(a," enr ^)[24) » Y(a,i,k,el,af,e3,tp) 


D Σ 


el e2 tp€ DSetX(u,a,t,k,el,e2,af,tp) 
tp + [CTtoE3(a,r)/24] =t 
Dist(e2,af) + 0 
CTtoE3(a,r) + CTtoK(a,r) 


( MOGReq(af,a) « GTime(a, "enr ")/24) » X(u,a,i,k,el,e2,af.tp) 


22.222 


el e2 tpecDsSetY(a,ik,el,e2,af,tp) 
tp + [RCTtoE3(a,r)/24) =1 
Distíaf k) = O 
RCTtoE3(a,r) + RCTtol(a,r) 


(MOGReq(af,a) « GTime(a,"enr”) /24) « Y(a,i,k,e1,e2,af,tp) 


2,2, δ. δ. δ, 2. 


i εἰ εὖ e3 tpeDSetX(u,a,i af.el,e2,e3,1p) 
tp + [CTtoK(a,r)/24] =t 


( MOGReq(af.a) * GTime(a," offld^)|24) » X(u,a,i,af,e1,e2,e3,tp) 


< MOGEff * MOGCap(af,"nb”) V af,t. 


4) 


The conditional checks differ due to the multiple roles (origin, enroute 
and destination) that an airfield may serve. When the airfield is serving as an origin, the 
conditional check will ensure that only permissible airlift missions are allowed; i.e., valid 
combination of unit u, aircraft a, route r, and mission start time Ι. 

When the airfield is serving as an enroute, four conditional checks are 
required. The first conditional check ensures a valid combination of indices; i.e., 
permissible unit u, aircraft a, mission start time ¢ combination. The second conditional 
check ensures that the time of arrival at the airfield is the time period of concern. The 
third conditional check ensures that the aircraft come from a different airfield and serves 
the purpose of preventing double counting. The fourth conditional check ensures that the 
time taken to reach the airfield of concern is not equal to the time required to reach the 
destination (or origin) airfield. This last check also ensures that proper aircraft ground 
time is used in the MOG computation. For example, if the time taken to reach enroute 
e3 and destination K are the same, this implies that the aircraft has reached its destination 
and no enroute airfield capacity should be consumed. The corresponding aircraft off-load 
MOG consumption is computed when the airfield is considered as a destination node. 

When serving as a destination airfield, two conditional checks are 
required. The first conditional check ensures a valid combination of indices. The second 
conditional check ensures that the time of arrival at the airfield is the time period of 


concern. 
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V. PERFORMANCE, ANALYTIC INSIGHTS AND LIMITATIONS 


A. PERFORMANCE 

The performance of the optimization model on the IBM RS6000 model 590 Unix 
workstation is relatively fast. For a sample problem with 20 units, 7 aircraft types, 17 
airfields and 30 time periods, GAMS/XA and GAMS/OSL generated and solved the LP 
problem in about 90 seconds. There were 8273 permissible decision variables and 6349 
constraint equations after all the problem reduction using dynamic sets were performed. 
The optimization model takes about 10 times longer when the same problem is solved 
on a personal computer (486/33 AT machine with 32MB of RAM). The total generation 
and solve times on the personal computer is about 15 minutes. 

The manner in which the optimization model is implemented reduces the amount 
of external data processing and data entry time (as most data can be entered in the format 
in which they are made available). For example, with readily available airfields' 
geographical location and aircraft's cruise speed data, the model computes the cumulative 
time a particular aircraft takes to reach a designated airfield along a route. The data 
entry time for the above sample problem is about one and a half hours. 

Thus, it can be seen that the overall response time of the model (under 2 hours) is 


relatively fast; this is a strength of the model. 
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B. ANALYTIC INSIGHTS 
Broad insights into the strategic airlift system can be gained from the optimization 
model. For any particular set of input data, the insights gained include answers to the 


following mobility questions: 


€ Adequacy of Assets: Indicate whether aircraft and airfield assets are adequate for 
a deployment scenario. 


@ Impact of Any Shortfall in Airlift Capability: Allow the analyst to assess the 
impact by reporting the number of troops and tons of equipment of each unit that 
are airlifted late or could not be airlifted. 


@ System "bottle-necks": Indicate which part (aircraft or airfields) of the airlift 
system is restricting the flow of troops and equipment. These reports can also give 
ideas on how to improve the airlift system output. For example, increasing the 
material handling equipment in those airfields that support a large number of 
aircraft (high MOG consumption) can increase aircraft handling capacity. An 
analyst can also gauge the optimal number of CRAF aircraft to mobilize. For 
example, if the system is airfield constrained then mobilizing additional CRAF 
assets may not help the system; on the contrary, it may increase congestion and 
cost. 


@ Balance aircraft and airfields assets: This factor is important as there are limited 
military budgets and the USAF has to make a trade-off between buying aircraft and 
keeping airfields open. Ideally, the combination of assets should be such that 
aircraft utilization rates and airfields MOG consumption for most deployment 
scenarios are close to their published capacities. This is to ensure a good aircraft 
to airfield ratio. 


e Fleet Mix Studies: By running different deployment scenarios, an analyst is able 
to assess which assets-mix option is the best overall option for supporting the 
different regional contingency plans. 


C. MODEL LIMITATIONS 

Two limitations of the optimization model are the inability to handle local traffic 
congestion and low fidelity. The true impact of these limitations are not quantifiable at 
the moment, and further validation and research are needed before it can be ascertained. 
These two limitations are caused by a mismatch between the time resolution (in days) and 
aircraft flying time and ground time (in hours). A more refined time resolution was not 
used in this research as initial investigation indicated that the problem size is large and 
requires data aggregation. A more refined time resolution though desired (more time 
periods for the same number of days) may increase the LP size beyond GAMS or the 


solvers' capability. This time resolution issue, however, warrants further investigation. 


1. Local Traffic Congestion 
Although the aircraft handling capacity of each airfield is observed by the 
model, airlift missions may still be routed in a manner that causes local congestion. For 
example, all aircraft routed through an airfield on a particular day could arrive within 
a small time window instead of being spread over the whole day. In reality, this would 
cause local congestion, even though the model's representation of aircraft handling 
Capacity is observed. 
2. Low Fidelity 
The arrivals of aircraft at the origin and destination airfields are rounded-off 
to the nearest day (resolution of the time period). This round-off affects the fidelity and 


accuracy of the model. For example, rounding-off may increase the MOG consumption 
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of an airfield on one day and decrease the MOG consumption of the same airfield on the 
next day. Round-off results in lumpy flows through airfields, and hence, errs on the side 
of pessimism when MOG represents a binding constraint on the system flow. This is to 
be contrasted with the more optimistic smoothing fractional flow approach utilized by 


RAND in CONOP [Ref. 4]. 
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————  —— —OMT F -  . — 


VL CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The USAF strategic airlift system has been successfully modelled (although not in 
its entirety as use of tanker aircraft are not incorporated) as a multi-commodity, multi- 
period network flow model with a large number of side constraints. 

Broad insights (some of which are currently not available from existing models 
such as Thruput or MOM) into the strategic airlift system can be gained by using the 
model. Some of the insights that can be gained include identification of system "bottle- 
necks" and the impact of any shortfall in airlift capability. 

In addition to providing insights into the airlift system, the model also has relatively 
fast turnaround time. This virtue makes the model a useful tool 1n situations (e.g., initial 
planning phase of a major regional conflict development) where time ts of the essence 
and where quick answers are desired. Its ability to furnish quick answers without tying 
down a massive amount of manpower or computer resources is important in such time 


critical circumstances. 


B. RECOMMENDATIONS 


The following are some recommendations for model enhancement. 
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1. Incorporate Use of Tanker Aircraft 

One enhancement to the model is to incorporate the use of tanker aircraft in 
a deployment scenario. Such an enhancement would allow the airlift system to be 
examined from a much wider perspective. A possible way to incorporate the use of 
tanker aircraft is the "airbase in the sky" concept employed by RAND’s CONOP 
optimization model. This concept basically treats the aerial refuelling areas as enroute 
airfields and can be implemented in the existing model without major changes to GAMS 
code. A more realistic representation of aerial refueling would also accounts for 
"diverts", i.e., planned aerial refueling which are not executed for unplanned reasons. 


Incorporating this enhancement would require major research and development. 


2. Increase MOG Fidelity 
For a more accurate modelling approach, it is essential that "ground services 
MOG" be distinguished from "parking space MOG". This distinction is necessary as 
aircraft inventoried at the origin and destination airfields require only "parking space 
MOG" and not "ground services MOG" once they have loaded or unloaded their cargo. 
This point has been covered in the Assumptions section in Chapter II. If data were 
available, it would be easy to formulate MOG constraints for each of several ground 


services. 


3. Investigate the Feasibility of Improving Time Resolution 
Another possible enhancement to the model is to improve the resolution of 


the time period (e.g., tume could be discretized into 12 hour or smaller blocks instead of 
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days). Such an improvement if feasible would certainly improve the fidelity of the model 
and alleviate the two limitations of the models cited in Chapter V. However, improving 
the resolution of the time period does has a high price in terms of increasing the size of 
the LP. It is therefore prudent to assess the relative merits before embarking on such a 


move. 


4. Validate MOG Efficiency Factor 
The MOG efficiency factor is introduced as a discount factor to offset the 
effect of random aircraft ground times (which implies that an airfield cannot be planned 
to be fully utilized). The appropriate efficiency factor however, is presently unknown 
and needs to be investigated. One possible approach is to employ a simulation model to 


help obtain a realistic efficiency factor. 


5. Stochastic Modelling 
Aircraft reliability and ground times at on-load, off-load, and enroute airfields 
are inherently random aspects of a strategic airlift system. These stochastic factors can 
significantly affect the performance of an airlift system, particularly when infrastructure 
(e.g., airfield capacity) represents a binding constraint. One possible enhancement to the 
deterministic model is to apply stochastic optimization to develop a strategic airlift model 
in which aircraft reliability and ground times are modeled as random variables with 


known distribution. 
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APPENDIX A. LISTING OF DATA SOURCES 


This appendix is a listing of the data sources. The layout is in the order of the 
entity name (as in the GAMS formulation) followed by a short description and the source 
from which the data originate. 


A. MOVEMENT REQUIREMENT DATA 


Name 


MoveUE(u,i,k) 
CargoP(u,cc) 
MovePAX(u,i,k) 
ALD(u) 

RDD(u) 


B. PENALTY DATA 


Name 


LatePen(u,c) 
NoGoPen(u,c) 
MaxLate 


C. CARGO DATA 


Name 


UESqFt(u) 
PAXWIt 


D. AIRCRAFT DATA 


Name 


Supply(a,t) 
ACSize(a, ab) 
ACCargo(a,cc) 
MaxPAX(a) 
PAXSqft(a) 
ACSqFt(a) 


Description 


Movement requirement for equipment 
Proportion of cargo classes (bulk, over, out) 
Movement requirement for troops 

Available load date 

Required delivery date 


Description 


Lateness penalty 
Non-delivery penalty 
Maximum allowed lateness for delivery 


Description 


Cargo space taken up by equipment 
Average troop weight 


Description 


New aircraft made available 
Aircraft body type 

Aircraft-cargo compatibility 
Aircraft troop carriage capacity 
Cargo space consumed by a soldier 
Cargo space available 
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Source 


Scenario 
scenario 
scenario 
Scenario 
Scenario 


Source 


analyst's input 
analyst's input 
analyst's input 


Source 


scenario 
AMC Load Planning Guide 


AFR 76-2, AMCP 55-41 
aircraft-9, 55-41 
AMC Load Planning Guide 


D. AIRCRAFT DATA (continued) 





Name Description 

LoadEff(a) Cargo space loading efficiency 
Spd(a) Aircraft block speed 
GTime(a,ga) Aircraft ground time 

URate(a) Established utilization rate 


ACRange(a,rangeband) Aircraft range points 
ACLoad(a,rangeband) Aircraft load at range point 


E. AIRFIELD DATA 


Name Description 

Crd(af, coords) Decimal airfield coordinates 
MOGCap(af, ab) MOG capacity by aircraft type 
MOGEff MOG efficiency factor 


F, AIRCRAFT-ROUTE DATA (Dynamic Sets) 


Name Description 


VRouteX(a,i,k,el,e2,e3) Aircraft-delivery route compatibility dynamic set 
VRouteY(a,i,k,el,e2,e3) Aircraft-recovery route compatibility dynamic set 


Source 


analyst's input 
AMC 

AMC 

AMC 
aircraft-9, 55-41 
aircraft-9, 55-41 


Source 


JFAST GEOFILE database 
AMC 
analyst's input 


Source 


analyst's input 
analyst's input 


The above are the original data that are required by the program. All other data 


are computed by the GAMS model itself. 


GAMS 2.325.064 


APPENDIX B. GAMS FORMULATION LISTING 


09/07/94 09:58:22 
Strategic Airlift Asset Optimization Model Spring 1994 


AIX RS/6000P 


2 * Thesis by: Cpt Lim, Teo-Weng 
3 * Advisors: Professor Richard E. Rosenthal 
4 * Dr. David P. Morton 
5 * Professor Gerald G. Brown 
6 * Naval PostGraduate School 
7 * 
8 *--------------- GAMS AND DOLLAR CONTROL OPTIONS --------------- 
9 
14 OPTIONS 
15 LIMCOL = 0, LIMROW = O , SOLPRINT = OFF, DECIMALS = 2 
16 RESLIM = 18000, ITERLIM = 400000, OPTCR = 0.1 , SEED = 3141; 
18 
19 αφ INDICES DECLARATION +++++++++++++++++++++ 
20 
21 SEIS 
22 
23 U units 
INCLUDE /home/limt/unitname.dat 
25 / 
26 UnitA USMC 3 MAB 
27  UnitB USMC 1 MAF 
28 UnitC USMC 2 MAB 
29 1 UnitD USA 4 Light Inf 
30 UnitE USA 2 Mech 
31 UnitF USA 25 AAD 
32 UnitG 24 AlOs 
33 UnitH 24 Flés and 24 F15s 
34  UnitI 24 AlOS 
35  UnitJ 24 F16s and 24 F15s 
36 UnitK 24 AlOs 
375 Unit D 24 F16s 
38 UnitM 24 F15s 
39 UnitN 24 AlOs 
40 UnitO 24 F16s 
41 UnitP 24 F15s 
42 UnitQ 24 F16s and 24 F15s 
43  UnitR 24 Fl6s and 24 Ε155 
44 Units 24 F16s and 24 F15s 
45  UnitT 24 Fl6s and 24 F15s 
46 
47 
48 A aircraft 
INCLUDE /home/limt/acname.dat 
50 
51 C5 
52 ΕΠ 
53 C141 
54 C130 
55 747P 
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56 747C 
57 DC10 
58 / 
59 
60 CRAF (A) cividjyan aircraft 
INCLUDE /home/limt/craf.dat 
em / 
63 747P 
64 747C 
65  DC10 
ee / 
67 
68 MILAC (A) military aircraft 
INCLUDE /home/limt/milac.dat 
70 / 
gu ον 
725 C17 
73 C141 
74 C130 
75 / 
76 
77 AB aircraft body type 
78 * Aircraft are classified as WB (widebody), NB (narrow body), 
79 * GM (ground manuverable), or TT (tactical). 
80 / WB, NB, GM, TT / 
81 
82 GA ground activities 
83 * The different types of ground activities an aircraft goes 
84 * through are on-load (ONLD), enroute turnaround (ENR) and 
85 * off-load (OFFLD). 
86 / ONLD, ENR, OFFLD / 
87 
88 T time periods in days 
INCLUDE /home/limt/periods.dat 
90 * Time period should cover up till the last RDD plus 
91 * the maximum allowed lateness 
92 / 
93 DAY1 * DAY30 
94 
95 
97  ALIAS(T,TP); 
98 
99 SETS 
100 
15 C cargo type / UE, PAX / 
102 * UE stands for unit equipment and PAX for troops 
103 
104 τε cargo class / BULK, OVER, OUT / 
105 
106 AF airfields in general 
INCLUDE /home/limt/afname.dat 
108 / 
109 QFQE Mildenhall 
110 TMKH Pope AFB 
111 TYFR Ramstein AB 
2802. XDAT Travis AFB 
113 FFTJ Dhahran 
114  UGZX Riyadh 
115  PKVV Masirah 
116 KNMD Hickam 
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117  FXSB Elmendorf AFB 

118 WWYK Tinker AFB 

119 NNGX Boston Logan 

120 UTKY San Francisco Int 
121 FJXT Dover AFB 

122 ZNRE Yokota 

123 HTDS London Gatwick 

124  FGDC Diego Garcia NAF 

125  XBGX Torrejon 

126 / 

127 

129 I(AF) origin airfields 
INCLUDE /home/limt/orgname.dat 
131 / 

132 QFQE Mildenhall 

133 TMKH Pope AFB 

134  TYFR Ramstein AB 

135  XDAT Travis AFB 

136 / 

137 

138 K (AF) destination airfields in theater 
INCLUDE /home/limt/destname.dat 
140 

141  FFTJ Dhahran 

142 UGZX Riyadh 

143 PKVV Masirah 

144 / 

145 

146 El (AF) first set of enroute airfields 
INCLUDE /nome/limt/elname.dat 
148 

149  QFQE Mildenhall 

150 TMKH Pope AFB 

151  TYFR Ramstein AB 

152  XDAT Travis AFB 

153 KNMD Hickam 

154 FXSB Elmendorf AFB 

155 WWYK Tinker AFB 

156 NNGX Boston Logan 

157 UTKY San Francisco Int 
158 ΕΟΧΤ Dover AFB 

159 / 

160 

161 E2 (AF) second set of enroute airfields 
INCLUDE /home/limt/e2name.dat 
163 / 

164 QFQE Mildenhall 

165 TYFR Ramstein AB 

166 FFTJ Dhahran 

167 UGZX Riyadh 

168 PKVV Masirah 

169 KNMD Hickam 

170 ZNRE Yokota 

171  HTDS London Gatwick 

172 f 

173 

174 E3 (AF) third set of enroute airfields 
INCLUDE /home/limt/e3name.dat 
176 / 

177 TYFR Ramstein AB 
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178 FFTJ Dhahran 

179  UGZX Riyadh 

180 PKVV Misarah 

181  FGDC Diego Garcia NAF 

182 ΧΒΟΧ Torrejon 

183 / 

184 ; 

185 

186 SET VROUTEX(A,I,K,E1,E2,E3) indicate aircraft-delivery route 

τ compatibility ; 

187 * Used to control aircraft-route matching; i.e., allow military 
188 * aircraft to fly on military routes and civilian aircraft to fly 
189 * on civilian routes. Members of SET indicate compatibility. 
INCLUDE /home/limt/vroutex.set 

191 * Member indicates aircraft-route compatibility 

192 * Substitute destination name for El, E2, E3 when such enroute 
193 * airfields are not needed. 

194 VRouteX(MILAC,"TMKH","FFTJ","WWYK","FFTJ","FFTJ") = yes; 

195 VRouteX (MILAC, "TMKH", "UGZX", "WWYK", "UGZX", "UGZX") = yes; 

196 VRouteX (MILAC, "TMKH","FFTJ", "TYFR","FFTJ","FFTJ") = yes; 

197 VRouteX (MILAC, "TMKH", "UGZX", "TYFR", "UGZX", "UGZX") = yes; 

198 VRouteX (MILAC, "TMKH", "FFTJ","QFQE","FFTJ","FFTJ") = yes; 

199 VRouteX (MILAC, "TMKH", "UGZX", "QFQE", "UGZX", "UGZX") = yes; 

200 VRouteX(CRAF,"TMKH","FFTJ","TMKH","TYFR","FFTJ") = yes; 

201 VRouteX (CRAF, "TMKH", "UGZX","TMKH", "TYFR", "UGZX") = yes; 

202 VRouteX(CRAF, "TMKH","PKVV","TMKH", "TYFR","PKVV") = yes; 

203 VRouteX (CRAF, "TMKH", "FFTJ", “TMKH", "HTDS","FFTJ") = yes; 

204 VRouteX(CRAF, "TMKH", "UGZX","TMKH", "HTDS","UGZX") = yes; 

205 "VRouteX(CRAF, "TMKH", "PKVV", "TMKH", "HTDS", "PKVV") = yes; 

206 VRouteX (MILAC, "XDAT","FFTJ", "WWYK","FFTJ","FFTJ") = yes; 

207 VRoutexX (MILAC, "TMKH", "FFTJ", "WWYK", "FFTJ","FFTJ") = yes; 

208 VRouteX (MILAC, "XDAT", "UGZX", "WWYK", "UGZX", "UGZX") = yes; 

209 VRoutexX (MILAC, "TMKH", "UGZX", "WWYK", "UGZX", "UGZX") = yes; 

210 VRouteX(MILAC,"XDAT","FFTJ","TYFR","FFTJ","FFTJ") = yes; 

211 VRouteX (MILAC, "TMKH", "FFTJ", "TYFR", "FFTJ","FFTJ") = yes; 

212 VRouteX(MILAC, "XDAT","UGZX","TYFR","UGZX","UGZX") = yes; 

213 VRouteX (MILAC, "TMKH", "UGZX","TYER","UGZX","UGZX") s yes; 

214 VRouteX(MILAC, "XDAT","FFTJ","QFQE","FFTJ","FFTJ") = yes; 

215 VRouteX(MILAC, "TMKH","FFTJ","QFQE","FFTJ","FFTJ") = yes; 

216 VRouteX(MILAC, "XDAT","UGZX","QFQE","UGZX","UGZX") = yes; 

217 VRouteX (MILAC, "TMKH", "UGZX", "QFQE", "UGZX", "UGZX") = yes; 

218 VRouteX (CRAF, "XDAT", "FFTJ", "NNGX", "TYFR","FFTJ") = yes; 

219 VRouteX(CRAF, "XDAT","UGZX","NNGX","TYFR","UGZX") = yes; 

220 VRouteX(CRAF, "XDAT", "PKVV","NNGX", "TYFR","PKVV") = yes; 

221  VRouteX(CRAF,"XDAT","FFTJ","NNGX","HTDS","FFTJ") = yes; 

222 VRouteX(CRAF, "XDAT", "UGZX", "NNGX", "HTDS", "UGZX") = yes; 

223  VRouteX(CRAF,"XDAT","PKVV","NNGX","HTDS","PKVV") = yes; 

224 SURGUteX(A, TYFR","FFTJ","TYFR","TYFR","FFTJ") = yes; 

225 VRoutex(A, "TYFR", "UGZX"."TYFR", "TYFR", "UGZX") = yes; 

226 VRoutex(A, "TYFR","PKVV", "TYFR", "TYFR","PKVV") = yes; 

227 VROUteX(A, "QFQE","FFTJ", "TYFR","FFTIJ","FFTJ") = yes; 

228 VRouteX(A, "QFQE", "UGZX","TYFR", "UGZX","UGZX") = yes; 

229  VRouteX(A, "QFQE","PKVV", "TYFR","PKVV","PKVV") = yes; 

230 VRouteX(CRAF,"XDAT",K,"FXSB","ZNRE","FGDC") = yes; 

231 VRouteX(CRAF, "XDAT",K,"XDAT","ZNRE","FGDC") = yes; 

232  VRouteX(MILAC,"TMKH","PKVV","WWYK","UGZX","UGZX") - yes; 

233 VRouteX(MILAC, "TMKH","PKVV","QFQE","UGZX","UGZX") = yes; 

234 VRouteX(MILAC, "TMKH", "PKVV","TYFR","UGZX","UGZX") = yes 

235 

236 SET VROUTEY(A,I,K,E1,E2,E3) indicate aircraft-recovery route 
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* compatibility ; 


INCLUDE /home/limt/vroutey.set 

238 * Member indicates aircraft-route compatibility 

239 * Substitute destination name for El, E2, E3 when such enroute 
240 * airfields are not needed 

241 VRouteY(MILAC, "TMKH","FFTJ","WWYK","FFTJ","FFTJ") = yes; 
242 VRouteY(MILAC, "TMKH", "UGZX", "WWYK", "UGZX","UGZX") = yes; 
243 VRouteY(MILAC, "TMKH","FFTJ","TYFR","FFTJ","FFTJ") = yes; 
244 VRouteY(MILAC, "TMKH", "UGZX", "TYFR", "UGZX", "UGZX") = yes; 
245 VRouteY(MILAC, "TMKH","FFTJ", "QFQE","FFTJ","FFTJ") = yes; 
246 VRouteY (MILAC, "TMKH", "UGZX", "QFQE", "UGZX","UGZX") = yes; 
247 VRouteY(CRAF, "TMKH", "FFTJ","NNGX","FFTJ","FFTJ") = yes; 
248 VRouteY (CRAF, "TMKH", "UGZX", "NNGX", "UGZX","UGZX") = yes; 
249 VRouteY(CRAF,"TMKH","FFTI","WWYK", "FFTJ","FFTJ") = yes; 
250 VRouteY(CRAF,"TMKH", "UGZX", "WWYK", "UGZX","UGZX") = yes; 
251 "VRouteY(CRAF,"TMKH","FFTJ","TYFR","FFTJ","FFTJ") = yes; 
252 VRouteY (CRAF, "TMKH", "UGZX", "TYFR", "UGZX", "UGZX") = yes; 
253 VRouteY(CRAF, "TMKH", "FFTJ","QFQE","FFTJ","FFTJ") = yes; 
254 VRouteY(CRAF, "TMKH", "UGZX", "QFQE", "UGZX", "UGZX") = yes; 
255 VRouteY (CRAF, "XDAT",K, "FXSB","ZNRE","FGDC") = yes; 

256 VRouteY(CRAF,"XDAT",K,"KNMD","ZNRE","FGDC") = yes; 

257 VRouteY(CRAF, "XDAT",K, "UTKY","ZNRE","FGDC") = yes; 

258 VRouteY (MILAC, "TMKH", "PKVV", "WWYK", "UGZX","PKVV") = yes; 
259 VRouteY (MILAC, "TMKH","PKVV","QFQE", "UGZX","PKVV") = yes; 
260 VRouteY(MILAC,"TMKH","PKVV","TYFR","UGZX","PKVV") = yes 
261 

262 SET COORDS coordinates / LAT, LON /; 

263 

264 *+++++++++++++++++ MOVEMENT REQUIREMENT DATA +++++++++++++++++++ 
265 

266 PARAMETER MOVEUE (U, I,K) units' equipment movement requirements 
INCLUDE /home/limt/moveue.dat 

268 * Entries are in terms of 100 tons and in the order of 
269 * unit name, origin airfield and destination airfield 
29079 

271 UnitA.XDAT.FFTJ = 186.77 

272 UnitB.TYFR.FFTJ = 6.92 

273 UnitC.XDAT.FFTJ w 63.23 

274  UnitD.TMKH.FFTJ = 114.94 

275 UnitE.TYFR.UGZX = 663.69 

276 UnitF.TMKH.UGZX = 238.27 

277 UnitG.TYFR.UGZX = 5.03 

278 UnitH.XDAT.UGZX = 10.13 

279 UnitI.QFQE.UGZX = 5.03 

280 UnitJ.XDAT.UGZX = 10.13 

281 UnitK.TYFR.UGZX = 5.03 

282 UnitL.TYFR.UGZX = 4.91 

283 UnitM.TYFR.UGZX = 5.22 

284 UnitN.TYFR.UGZX = 5.03 

285 UnitO.TYFR.PKVV = 4.91 

286 UnitP.TYFR.PKVV = 5.22 

287 UnitQ.TMKH.PKVV = 10.13 

288 UnitR.TMKH.UGZX = 10.13 

289 UnitS.TMKH.UGZX = 10.13 

290 UnitT.XDAT.UGZX = 10.13 

291 / 

292 : 

2. 

294 TABLE CARGOP(U,CC) 

295 *  Proportion of cargo types (bulk, over, out) belonging to 
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296. * unit u moving from i to k. Proportions for the different 


gee * types of cargo must add to 1. 
INCLUDE /home/limt/cargop.dat 
299 BULK OVER OUT 
300 UnitA 0.383 0.520 0.097 
301 UnitB 0.470 0.530 0.0 
302 Unitc 0.331 0.651 0.018 
303 UnitD 0.190 0.797 0.013 
304 UnitE 0.023 0.575 0.402 
305 UnitF 0.064 0.838 0.098 
306 UnitG 0.616 0.384 0.0 
307 UnitH 0.661 0.339 0.0 
308 UnitI 0.616 0.384 0.0 
309 UnitJ 0.661 0.339 0.0 
310 UnitK 0.616 0.384 0.0 
ag σας, 0.658 0.342 0.0 
312 UnitM 0.665 0.335 0.0 
313  UnitN 0.616 0.384 0.0 
314  UnitO 0.616 0.384 0.0 
315 UnitP 0.658 0.342 0.0 
316 UnitQ 0.661 0.339 0.0 
317 UnitR 0.661 0.339 0.0 
318 Units 0.661 0.339 0.0 
319 UnitT 0.661 0.339 0.0 
320 ; 

321 

322 PARAMETER MOVEPAX(U,I,K) units' troop movement requirements 


INCLUDE /home/limt/movepax.dat 
324 * Entries are in terms of 100 men and in the order of 
325 , unit name, origin airfield and destination airfield 
326 
327 UnitA.XDAT.FFTJ = 126.98 
328 UnitB.TYFR.FFTJ = 13.84 
329 UnitC.XDAT.FFTJ = 52.59 
330 UnitD.TMKH.FFTJ = 97.76 
331 UnitE.TYFR.UGZX = 199.37 
332 UnitF.TMKH.UGZX = 190.82 
333 UnitG.TYFR.UGZX = 3.92 
334 UnitH.XDAT.UGZX = 8.46 
335 UnitI.QFQE.UGZX = 3.92 
336 Unitd.XDAT.UGZX = 8.46 
337 UnitK.TYFR.UGZX = 3.92 
338 UnitL.TYFR.UGZX = 4.74 
339 UnitM.TYFR.UGZX = 3.72 
340 UnitN.TYFR.UGZX = 3.92 
341 UnitO.TYFR.PKVV = 4.74 
342 UnitP.TYFR.PKVV = 3.72 
343 UnitQ.TMKH.PKVV = 8.46 
344 UnitR.TMKH.UGZX = 8.46 
345 UnitS.TMKH.UGZX = 8.46 
346 UnitT.XDAT.UGZX = 8.46 
347 
348 
349 PARAMETER ALD (U) available load date for unit u 

INCLUDE /home/limt/ald.dat 
351 
352  UnitA 8 
353 UnltB 16 
354 UnitC 3 
355 UnitD 21 
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PARAMETER RDD(U) 
/home/limt/rdd.dat 


SCALAR MAXLATE 
/home/limt/maxlate.dat 


P 


Ph bh P Q) O O) O Ñ N) Ñ. Ñ) eH Fe Q 


0 0 Ὁ) Ὁ —J Οἱ 6ἱ OA OA (n un Ma qe 


required delivery date for unit u 


maximum allowed lateness in days 


* Maximum allowed lateness for delivery in days 


TABLE LATEPEN (U,C) 
/home/limt/latepen.dat 


lateness penalty per day by unit and cargo 


* Late penalty for unit equipment should be measured in terms 


* of per 100 ton per day. 


Late penalty for PAX should be 


* measured in terms of per 100 PAX per day. 


356 UnitE 
357 UnitF 
358 UnitG 
359 UnitH 
360 UnitI 
361 UnitJ 
362 UnitK 
363 UnitL 
364  UnitM 
365 UnitN 
366 UnitO 
367 UnitP 
368 UnitQ 
369 UnitR 
370 UnitS 
371 BURIEtT 
37299 

373 

374 
INCLUDE 

376 / 

377 UnitA 
378  UnitB 
379 Unitc 
380 UnitD 
381  UnitE 
382 UnitF 
383  UnitG 
384  UnitH 
385  UnitI 
386 UnitJ 
387 UnitK 
388 UnitL 
389  UnitM 
390 UnitN 
391 UnitO 
392 UnitP 
393  UnitQ 
394  UnitR 
395 UnitS 
396 UnitT 
397 

398 

399 
INCLUDE 

401 

40205 1 
404 x: 

405 

406 
INCLUDE 

408 

409 

410 

411 

412 UnitA 
413 UnitB 
414 UnitC 
415 UnitD 
416 UnitE 


a oa 


.005 
.14 

.016 
.009 
. 002 


OOOOO 
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217 UnitF 0.004 0.005 
418 UnitG 0.20 0.255 
419 UnitH 0. 157 0.236 
420 UnitI 0.199 0.255 
421  UnitJ 0.197 0.236 
422 UnitK 0.199 0.255 
423 UnitL 0.204 0.217 
424  UnitM 0.192 0.269 
425  UnitN 0.199 0.255 
426 UnitO 0.204 0.211 
427 UnitP 0.192 0.269 
428 UnitQ 0.197 0.236 
429 UnitR 0.197 0.236 
430 UnitS 0.197 0.236 
431 UnitT 0.197 0.236 
432 ; 

433 


434 PARAMETER NOGOPEN(U,C) penalty for not fulfilling requirement ; 
435 *  NOGOPEN should be larger than (MAXLATE x LATEPEN) for the 
436 * respective unit and cargo type. This is for consistency or 
437 * else the LP may choose not to send instead of sending the 
438 * cargo late. 


INCLUDE /home/limt/nogopen.dat 
439 NOGOPEN(U,C) = 2 * MAXLATE * LATEPEN(U,C) 
440 ; 


441 Φε. ΡΕ CARGO DATA +++++++++++++++++++++++++ 
442 * 
443 PARAMETER UESQFT (U) ave. cargo space (in 1000 sq. ft.) per 100 


444 * tons of unit u's cargo 
INCLUDE /home/limt/uesqft.dat 
446 * Entries are in terms of 1000 sq. ft. per 100 ton 
447 / 

448  UnitA 47 

449  UnitB 2:6 

450 UnitC 2 9 

451 UnitD 3.0 

452 UnitE 2.4 

453 UnitF 3. 1 

454  UnitG 2.4 

455 UnitH 2.4 

456 UnitI 2.4 

457 UnitJ 2.4 

458 UnitK 2.4 

459 UnitL 2.4 

460 UnitM 2.4 

461 UnitN 2.4 

462 UnitO 2.4 

463 UnitP 2.4 

464 UnitQ 2.4 

465 UnitR 2.4 

466 UnitS 2.4 

467 UnitT 2.4 

468 

469 


470 SCALAR PAXWT ave. wt (in 100 tons) of 100 troops / 0.2 /; 

471 * Average PAX weight is set at 400 lbs or 0.2 stons inclusive of 
472 * personal equipment. Therefore, average weight of 100 PAX is 
473 * 0.2 hundred stons. 

474 

475 *++++++++++++++++++++++ AIRCRAFT DATA ++++++++++++++++++++++++++ 
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476 
477 


PARAMETER SUPPLY (A,T) no. of new aircraft made available on day T 


INCLUDE /home/limt/acsupply.dat 


479 * SUPPLY(A,T) represents additional number of aircraft (A) that 
480 * are first made available for use at time period (T). 
481 / 

492 C5.DAY1 = 10 

483 C5.DAY6 = 30 

484  C5.DAY11 - 40 

485 C17.DAY1 = 2 

486 C17.DAY6 = 6 

487  C17.DAY11 = 8 

488  C141.DAY1 = 20 

489 C141.DAY6 = 60 

490 C141.DAY11 = 80 

491 C130.DAY1 = 20 

492 C130.DAY6 = 60 

493 C130.DAY11 = 80 

494 747P.DAY6 = 10 

495 747P.DAY11 = 30 

496 747P.DAY16 = 20 

497 747C.DAY6 = 5 

498 747C.DAY11 = 25 

499 747C.DAY16 = 50 

500 DC10.DAY6 = 5 

501 DC10.DAY11 fe y25 

502 DC10.DAY16 = 50 

503 / 

504 : 

505 

506 SET ACSIZE(A, AB) aircraft body type 
INCLUDE /home/limt/acsize.dat 

508 / 

509. C5.WB 

S10 οι GM 

511 C141.NB 

512 C130. TI 

513 747P.WB 

514 747C.WB 

515 DCIO.WB 

516 / 

517 Š 

518 

519 TABLE ACCARGO((A, CC) aircraft and cargo class matching 
INCLUDE /home/limt/accargo.dat 

521 * An entry of 1.0 represents compatibility 
522 BULK OVER OUT 

5231008 1.0 1.0 1.0 

524 C17 1:50 ο 1.0 

525 C141 10 1.0 0.0 

526. Cl30 το que ο. ο 

527  747P 0.0 0.0 0.0 

528 747C 170 0.0 0.0 

5294 DCIO 0.0 0.0 0.0 

5305097 

531 

532 PARAMETER MAXPAX (A) troop carriage capacity (in 100s) of aircraft 
INCLUDE /home/limt/acmaxpax.dat 

534 * Entries are in terms of 100s of men 
53507 
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536 C5 0 73 

537 C17 0.54 

538 C141 1.43 

S39 CI30 0.44 

540 747P 4.02 

541 747C 0.0 

542 DC10 ο 

543 / 

544 ; 

545 

546 PARAMETER PAXSQFT (A) cargo space (in 1000 sq. ft.) occupied by 

X 100 PAX 

INCLUDE /home/limt/paxsqft.dat 

548 * Entries are in terms of 1000 sq. ft. per 100 PAX 

549 * C5 PAX passenger seats upstairs and does not consume cargo space 
550 * CRAF PAX haulers have 0 sq. ft. per PAX and have 0 sq. ft. of 
551 * cargo space. 

552 / 

553 C5 0-0 

554 C17 0.7 

555 C141 oc 

556 C130 0r 7 

557 747P 0.0 

558 747C 0.0 

559 DC10 0.0 

560 / 

561 ; 

562 

563 PARAMETER ACSQFT (A) alrcraft cargo space (in 1000 sq. ft.) 
INCLUDE /home/limt/acsqft.dat 

565 * Entries are in terms of 1000 sq. ft. 

566 * CRAF PAX haulers have O sq. ft. of cargo space. 

567 / 

568 C5 2:562 

569 C17 125593 

570 C141 1.043 

οτι οἱ 0.452 

572 747P 0.0 

573 747C 2.772 

574 DC10 0.0 

575 / 

576 E 

577 

578 PARAMETER LOADEFF (A) loading efficiency for aircraft type 
INCLUDE /home/limt/loadeff.dat 

580 * Loading Efficiency gives the proportion of cargo space that 
581 * could be actually filled. Loading Efficiency does not apply 
582 * to CRAF PAX haulers which have O sq. ft of cargo space. 
583 / 

584 C5 0.8 

585 C17 0.8 

586 C141 0.8 

587 C130 0.8 

588 747P 0.0 

589 747ς Je. 0 

590 DC10 O O 

591 / 

592 | 

593 

594 PARAMETER SPD(A) aircraft block speed in knots 
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INCLUDE /home/limt/acspeed.dat 
/ 


596 

597 C5 419 

598 C17 440 

599 C141 397 

600 C130 215 

601 747P 450 

602 747C 450 

603 DC10 450 

604 / 

60S ; 

606 

607 TABLE GTIME(A, GA) aircraft ground activity duration in hrs 
INCLUDE /home/limt/gtime.dat 

609 ONLD ENR OFFLD 

610 ο 4225 225 3-25 

611. C17 2225 2.25 2.25 

612 C141 Que 2.25 2025 

613 C130 1.5 1x5 1.9 

614 747P 3.0 2-0 3.0 

615 747C 4.0 2.0 3.0 

616 DC10 3.0 20 3.0 

GITE 

618 

619 PARAMETER URATE (A) aircraft planned utilization rate 

620 * Gives the planning norm for aircraft utilization. 
INCLUDE /home/limt/urate.dat 

627 * Entries are in 100 hrs per aircraft per day 

623 77 

624 C5 0.109 

625.7017 0-152 

626 C141 0.102 

627 C130 0.054 

628 747P 0.10 

629 747C 0.10 

630 DC10 0.10 

631 / 

632  ; 

633 

634  *--LAL AIRCRAFT RANGE-PAYLOAD DATA +++++++++++++++T++++ 
635 3 The following tables allow for the computation of critical 
636 * payload for a particular route. The entries in ACRANGE and 
637 *  ACLOAD (for a rangeband and aircraft pair) give the range and 
638 * respective critical payload data for the aircraft. These data 
639 * are used to interpolate the critical payload for a given range. 
640 

641 SET RANGEBAND / 1*7 /; 

642 

643 TABLE ACRANGE (A, RANGEBAND) aircraft range points 
INCLUDE /home/limt/acrange.dat 

645 1 2 3 4 5 6 7 

646 CS 2000 2500 3000 3500 4000 4500 10000 

647 C17 2400 3400 4200 10000 0 0 0 

648 C141 2000 2500 3000 3500 4000 4500 10000 

649 C130 2000 4200 4800 10000 0 0 0 

650 747P 4450 6200 10000 0 0 0 0 

651 747C 2425 3000 3500 6000 0 0 0 

652 DC10 3000 3500 6000 10000 0 0 0 

657. 

654 
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655 TABLE ACLOAD (A, RANGEBAND) 


aircraft payloads (in 100 tons) at range 


INCLUDE /home/limt/acload.dat 

657 * Entries are in terms of 100 tons 

658 T 2 3 4 5 6 7 
ος. 5 1 127 ο ος .0:873 0.757 0.630 0.506 0 
660 C17 0.800 0.650 0 0 0 0 0 
661 C141 0.344 0.299 0.259 0.203 0.161 0.114 0 
662 C130 0.215 0.085 0 0 0 0 0 
663 747P 0.64 0 0 0 0 0 0 
664 747C 1.286 -1.216 1.116 0 0 0 0 
665 DC10 0.44 0 0 0 0 0 0 
666 ; 

667 

668 PARAMETER SLOPE (A, RANGEBAND) slopes for interpolation; 
669 * Slopes give the change in payload from rangeband i-1 to 
670 * rangeband i. 

671 

672 SLOPE(A,RANGEBAND) = 

673 OS(ORD(RANGEBAND) EQ 1) + 

674 ( (ACLOAD (A, RANGEBAND- 1) - ACLOAD (A, RANGEBAND) ) / 

675 (15 (ACRANGE (A, RANGEBAND) EQ 0) + 

676 (ACRANGE (A, RANGEBAND) - ACRANGE (A, RANGEBAND- 1) ) 
677 $ (ACRANGE (A, RANGEBAND) GT 0)) )S(ORD(RANGEBAND) GT 1); 
678 

679 TW+++++++++++++++++++++++ AIRFIELD DATA +++++++++++++++++++++++++ 
680 

681 TABLE CRD(AF,COORDS) airfield coordinates in decimals 
INCLUDE /home/limt/afcoord.dat 

683 LAT LON 

684 QFQE 52.4 0 5 

685 TMKH 35.2 -79.0 

686 TYFR 49.4 7.6 

687  XDAT 38.3 -121.9 

688 FFTJ 26. 3 50:2 

689 UGZX 24.7 46.7 

690 PKW 2077 58:9 

691 KNMD 21.3 -157.9 

692 FXSB 61.25 -149 8 

693 WWYK 35.4 -97.4 

694 NNGX 42.36 -71.01 

τιν 37.62 =122.4 

696 FJXT 359 1 -75.5 

697  ZNRE 35.8 139.4 

698 HTDS 5175 0.7 

699  FGDC -7.3 72.4 

700 ΧΒΟΧ 40.49 -3.46 

του ; 

702 

703 TABLE MOGCAP(AF,AB) MOG capacity at airfield 

704 * MOG in essence represents the airfield's aircraft handling 
705 * capacity. MOG capacity depends on the aircraft body 
706 * type. MOG figures are based on a host of factors such as 
707 * ramp space, fuel availability, MHE, servicing, reception 
708 * capability and contention with collocated combat units. 
INCLUDE /home/limt/mogcap.dat 

710 WB NB GM TT 

711  QFQE 5 9 9 

712  TMKH 17 32 32 

713  TYFR 8 22 22 

714 XDAT 23 41 41 
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TUS ΕΕΣ 10 14 15 

716  UGZX 6 TO 13 

717  PKVV l 3 3 

718 KNMD 30 41 41 

719 FXSB 8 16 16 

720  WWYK 8 16 16 

721  NNGX 7 7 7 

722 "UTKY 7 T 7 

723  FJXT 2 14 14 

724  ZNRE 10 20 20 

725  HTDS 3 12 12 

726 FGDC 2 5 5 

727  XBGX 10 12 1 2 

προ. 

729 

730 * set all tactical MOG availabilities to NB MOG/.75 

731 

ma MOGCAP(AF,'TT') = MOGCAP(AF, 'NB')/.75; 

733 

734 SCALAR MOGEFF efficiency factor for MOG usage 
INCLUDE /home/limt/mogeff.dat 

736 * MOG efficiency factor is to cater for the fact that MOG capacity 
737 * cannot be achieved due to the stochastic nature of ground time. 
738 / 0.9 / 

740  ; 

741 

742 X++++++++++++++++++++++ ROUTE PARAMETERS +++++++++++++++++++++++ 
743 

744 * The following equations are used to compute the distance between 
745 * any two airfields. 

746 

747 ALIAS (AF,AFP); 

748 

749 PARAMETERS FX(AF,AFP) result of sin cos argument of acos fcn 
750 DIST(AF,AFP) distance between airfields; 

τοι 

752 SCALAR DTOR conversion factor from degrees to radians ; 

753 DTOR = 0.017453293; 

754 

755 FX(AF,AFP) = 

756 COS (DTOR*CRD (AF, ’ LAT’) ) *COS (DTOR*CRD (AFP, ’ LAT’) ) * 

157 COS (DTOR* (CRD(AF,'LON')-CRD(AFP,'LON'))) + 

758 SIN(DTOR*CRD(AF,'LAT'))*SIN(DTOR*CRD(AFP,'LAT')); 

759 

760 DIST(AF,AFP)S(ORD(AF) NE ORD(AFP) ) 

761 - 3437* ((ARCTAN (SQRT (1- SQR(FX(AF,AFP)))/ 

762 FX(AF,AFP)))$(FX(AF,AFP) GT 0) + 

763 (3.141592653+ARCTAN (SQRT (1-SQOR(FX(AF, AFP))) 

764 /FX(AF,AFP)))$(FX(AF,AFP) LT 0)); 

765 

766 PARAMETER MAXLEG(I,K,E1,E2,E3) longest leg of each route ; 

767 MAXLEG(I,K,E1,E2,E3) = 

768 MAX (DIST (I,E1) , DIST(E1,E2) ,DIST(E2,E3) ,DIST(E3, Kou 
769 

770 PARAMETER MAXLOAD(A,I,K,E1,E2,E3) maximum payload on longest leg ; 
771 * Gives maximum payload of an aircraft type on a particular route 
J72 MAXLOAD(A,I,K,E1,E2,E3) $ VROUTEX(A, 1,K,E1,E2,E3) = 

773 ACLOAD (A, '1')$(MAXLEG(1,K,E1,E2,E3) LE ACRANGE (A, '1')) + 
774 (ACLOAD (A, '1')+SLOPE (A, '2') * (ACRANGE (A, '1') -MAXLEG (1,K,E1,E2,E3))) 
775 $ (MAXLEG(1,K,E1,E2,E3) GE ACRANGE(A,'1') AND 
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MAXLEG(I,K,E1,E2,E3) LE ACRANGE(A,'2')) 4 
(ACLOAD(A,'2') *SLOPE(A,'3') * (ACRANGE(A,'2') -MAXLEG(I,K,E1,E2,E23))) 
$(MAXLEG(I,K,El1,E2,E3) GE ACRANGE(A,'2') AND 
MAXLEG (1,K,E1,E2,E3) LE ACRANGE(A,'3')) + 
(ACLOAD (A,'3') «SLOPE(A,'4' ) * (ACRANGE(A, '3' ) -MAXLEG(I,K,E1,E2,E3))) 
S(MAXLEG(I,K,E1,E2,E3) GE ACRANGE(A,'3') AND 
MAXLEG(I,K,E1,E2,E3) LE ACRANGE(A,'4')) + 
(ACLOAD(A,'4')«SLOPE(A,'5') * (ACRANGE (A, ' 4“ ) -MAXLEG(I,K,E1,E2,E3))) 
S (MAXLEG (I, K, E1, E2, E3) GE ACRANGE(A,'4') AND 
MAXLEG(I,K,E1,E2,E3) LE ACRANGE(A,'5')) + 
(ACLOAD(A,'5') «SLOPE(A,'6') * (ACRANGE(A,'5') -MAXLEG(I,K,E1,E2,E23))) 
S (MAXLEG(I,K,E1,E2,E3) GE ACRANGE(A,’S’) AND 
MAXLEG(I,K,E1,E2,E3) LE ACRANGE(A,'6')) 
(ACLOAD (A,'6') «SLOPE (A, ' 7! ) * (ACRANGE (A, '6') -MAXLEG (I, K,E1,E2, E3))) 
$ (MAXLEG (1,K,E1,E2,E3) GE ACRANGE (A, '6') AND 
MAXLEG(I,K,El,E2,E3) LE ACRANGE(A,'7')) ; 


PARAMETER CTTOE1(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to travel from origin airfield 
* (I) to enroute airfield (El) along the route specified inclusive 
"or ground time. 
GITOEI(A,ISK.EI,E2,E3) SVROUTEX(A,I,K,E1,E2,E3) Ξ 
GTIME(A,'ONLD') + (DIST(I,E1)/SPD(A)); 


PARAMETER CTTOE2(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to travel from origin airfield 
* to enroute airfield (E2) along the route specified inclusive of 
* ground time. 
GDIOBE2(A,I,K;EI,E2,E3) | $SVROUTEX(A,I,K,El,E2,E3) = 
CTTOEl(A,I,K,E1,E2,E3) + 
(GTIME (A, "'ENR') + (DIST(E1,E2)/SPD(A)))SDIST(E1,E2); 


PARAMETER CTTOE3(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to travel from origin airfield 
* to enroute airfield (E3) along the route specified. 
EDRTOES(A,I,K,El,E2,E3) SVROUTEX(A,I,K,E1,E2,E3) ^ 
CTTOE2 (A, 1,K,E1,E2,E3) + 
(GTIME (A, 'ENR') + (DIST(E2,E3)/SPD(A)))$DIST(E2,E3); 


PARAMETER CTTOK(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to travel from origin airfield 
* to destination airfield (K) along the route specified. 
GITOK(A.I-K.ELI,E2;E3) SVROUTEX(A,I,K,EI,E2,E3) = 
STTOE3 (A, I,K,El,E2,E3) + 
(GTIME (A, 'ENR') + (DIST(E3,K)/SPD(A)))SDIST(E3,K) ; 


PARAMETER FLTTIME (A, 1,K,E1,E2,E3) flight time (in 100 hrs) ; 
* Gives the total flying hours consumed by an aircraft along a 
* specified route. 
FLTTIME(A,I,K,E1,E2,E3) SVROUTEX(A,I,K,E1,E2,E3) = 
(DIST(I,E1) + DIST(E1,E2) + DIST(E2,E3) 
+ DIST(E3,K)) / ( SPD(A) * 100 ) ; 


PARAMETER RCTTOE3(A,AF,AF,AF,AF,AF) ; 

* Cumulative time (in hours) taken to recover from destination 

* airfield (K) to enroute airfield (E3) along the specified route 

* inclusive of ground time. Note: (a) the route sequence is now 

* from K to E3 to E2 to El to I. (b) the ordering of indices is 

* preserved as such, as recommended by GAMS. 
RCTTOE3 (A, L K, E1,E2,ËE3)  SVROUTEY(A, l,K,El,E2,E3) - 
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GTIME(A,'OFFLD') + (DIST(E3,K)/SPD(A)) ; 


PARAMETER RCTTOE2(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to recover from destination 
* airfield (K) to enroute airfield (E2) along the specified route 
* inclusive of ground time. 
RCTTOE2 (A, 1,K,E1,E2,E3) SVROUTEY (A, 1,K,E1,E2,E3) = 
RCTTOE3 (A, I1,K,E1,E2,E3) + 
(GTIME (A, "ENR') « (DIST(E2,E3) /SPD(A)) ) SDIST(E2,E3) ; 


PARAMETER RCTTOEl(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to recover from destination 
* airfield (K) to enroute airfield (El) along the specified route 
* inclusive of ground time. 
RCTTOEl(A,I,K,ELl,E2,E3) SVROUTEY(A,I,K,E1,E2,E3) z 
RCTTOE2(A,I,K,E1,E2,E3) + 
(GTIME(A,'ENR/) « (DIST(E1,E2)/SPD(A)))SDIST(E1,E2); 


PARAMETER RCTTOI(A,AF,AF,AF,AF,AF) ; 
* Cumulative time (in hours) taken to recover from destination 
* airfield (K) to origin airfield (I) along the specified route 
* inclusive of ground time. 
RCTTOI(A,1I,K,E1,E2,E3) SVROUTEY(A,I,K,E1,E2,E3) = 
RCTTOE1 (A, 1,K,E1,E2,E3) + 
(GTIME(A,'ENR') « (DIST(I,E1)/SPD(A)))SDIST(I,E1); 


PARAMETER MOGREQ (AF,4A) normalized (to NB) MOG requirement ; 
* Normalized MOG requirement for an aircraft type to that of a 
* narrow body aircraft for the airfield. 
MOGREQ(AF,A) - MOGCAP(AF,'NB') / 
( MOGCAP (AF, 'WB') SACSIZE(A,'WB') + 
MOGCAP (AF,'GM') SACSIZE(A,'GM') + 
MOGCAP(AF,'TT')S$ACSIZE(A,'TT') + 
MOGCAP (AF,'NB')SACSIZE(A,'NB') ); 


Y +++++++++++++++++++++++ LP VARIABLES ++++++++++++++++++++++++++ 
INTEGER VARIABLES 

X(U,A,AF,AF,AF,AF,AF,T) number of airlift missions(aircraft) for 
* unit (U), by aircraft (A), on the route 
* specified with start time (T). 

Y(A,AF,AF,AF,AF,AF,T) no. of aircraft (A) recovered to base (I) 
* on the route specified with start 
* time (T). 

ALLOT(A,I,T) new aircraft made available allotted to base (I). 


H (A,I,T) number of aircraft (A) saved for use at time T+1, 
* at origin base (I), an end-of-period inventory. 


HP (A,K, T) number of aircraft (A) that will start recovery ; 
* from destination (K) in the next period instead of 
x the current period. 
POSITIVE VARIABLES 


TONSUE(U,A,I,K,E1,E2,E3,T) unit U's equipment (in 100 tons) 
* airlifted by aircraft (A) on the 
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E specified route with start time (T). 
TPAX(U,A,I,K,EIl,E2,E3,T) number (in 100s) of unit U's troops 

E airlifted by aircraft (A), on the 

* specified route with start time (T). 
UENOGO (U, I, K) unit U's equipment (in 100 tons) not airlifted 

τ from the origin base (I) to destination (K) 

* in the theater. 


PAXNOGO(U,I,K) number (in 100s) of unit U's troops not airlifted; 
* from the origin base (I) to destination (K) 
* in the theater. 


FREE VARIABLE 
2 total penalty incurred from late deliveries and for not 
* fulfilling delivery requirement (undelivered cargoes). 


k++++++++++++++++++++ DYNAMIC SET DECLARATION ++++++++++++++++++ 


SET DSETX(U,A,AF,AF,AF,AF,AF,T) ; 
Used to control the allowable combination of unit, alrcraft, 
route, and start time for the following decision variables: 
X IU A LI K El, E2,E3,T), TONSUE(U,A,1,K,E1,E2,E3,T) and 
TPAX(U,A,I,K,E1,E2,E3,T). A valid combination implies that the 
following be satisfied: Start time for delivery must be after an 
unit's ALD and delivery date must be before the RDD plus the 
maximum allowed lateness. There is movement requirement. 
Maximum aircraft payload for the route must be greater than 25% 
its full load capability; this is for cost efficiency. Supply 
of aircraft type to date must be greater than O. 
DSEIX(U,VROUTEX(A,I,K;EL,E2,E3),T) = YES Sl 
(ORD(T) GE ALD(U)) AND 
(ORD(T) LE (RDD(U) - ROUND(CTTOK(VROUTEX)/24) + MAXLATE)) 
AND (1$MOVEUE(U,I,K) OR 1$MOVEPAX(U,I,K)) 
AND (1$(MAXLOAD(VROUTEX) GT (0.25 * ACLOAD(A,"1")) )) 
and SUM (TPS(ORD(TP) LE ORD(T)), SUPPLY(A,TP)) 
us 


+ + * + X f + f + b 


SET DSETY(A,AF,AF,AF,AF,AF,T) : 
* Used to control the allowable combination of aircraft, 
* route, and start time for the decision variable 
* Y(U,A,I,K,El1,E2,E3,T). An allowable combination is one 
* in which the aircraft-route is compatible and that the supply 
* of aircraft type to date is greater than 0. 
DSETY(VROUTEY(A,I,K,E1,E2,E3),T) = YES S( 
SUM (TPS(ORD(TP) LE ORD(T)), SUPPLY(A,TP)) 
nh 


kt+tt+++++++++++++++ LP EQUATION DEFINITIONS +++++++++++++++++++++ 
EQUATIONS 


OBJFUN define objective function 

REQMTUE (U, 1,K) movement requirement for UE 
OUTREQMT (U, I,K) out-sized cargo requirement 
OVERREQMT (U, I,K) over-sized cargo requirement 
REQMTPAX (U, I,K) movement requirement for troops 
ACBALI (A,I,T) aircraft balance equation at origins 
ACALLOT (A, T) aircraft allocation balance equation 
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ACBALK (A,K,T) aircraft balance equation at destinations 
ACWEIGHT(U,A,I,K,E1,E2,E3,T) aircraft payload limitation 
ACSPACE(U,A,I,K,E1,E2,E3,T) aircraft cargo space limitation 
ACPAX(U,A,I,K,E1,E2,E3,T) aircraft PAX carriage limitation 
ACURATE (A) utilization rate constraint 

MOGUTILITY(AF,T) handling capacity of airfields ; 


γα γ.γ OBJECTIVE FUNCTION +++++++++++++++++++++++ 


OBJFUN.. 

Z =E= SUM ( (U,A,I,K,EL,E2,E3,T) SDSETX(U,A,I,K,E1,E2,E3,T), 
(TONSUE (U,A,I,K,E1,E2,E3,T)SACSQFT(A) * LATEPEN(U,'UE') * 
MAX(0, ORD(T) -ROUND(CTTOK(A,I,K,E1,E2,E3)/24) -RDD(U)) } 4 
(TPAX(U,A,I,K,EL,E2,E3,T) * LATEPEN(U,'PAX') * 
MAX(0, ORD(T) +ROUND(CTTOK(A, 1,K,E1,E2,E3)/24)-RDD(U)) )) 

+ SUM ( (U,I,K) S(1SMOVEUE(U,I,K) OR 1$MOVEPAX(U, I,K)), 
(UENOGO (U, I, K) *NOGOPEN(U,'UE')) + 
(PAXNOGO (U, I, K) *NOGOPEN(U, ' PAX! )) }; 


f+++++++++++++++ DEMAND SATISFACTION CONSTRAINTS ++++++++++++++++ 


Total tonnage (in 100 tons) of unit equipment delivered to the 
theater by all aircraft types flying all routes over all time 
periods plus total tonnage (in 100 tons) of equipment not 
delivered must be greater than or equal to the movement 
requirement for unit equipment. 


+ f + ++ 


REQMTUE(U,I,K) $MOVEUE(U,I,K).. 
SUM ( (A,E1,E2,E3) S$ACCARGO(A,'BULK'), 
SUM ( T SDSETX(U,A,1I,K,E1,E2,E3,T), TONSUE(U,A,1,K,E1,E2, 5375 0 
+ UENOGO(U,I,K) =G= MOVEUE(U,I,K) ; 


Total tonnage (in 100 tons) of unit equipment carried by out- 
sized cargo capable aircraft plus tonnage of unit equipment 
not delivered (UENOGO) must be greater than or equal to the 
movement requirement for out-sized cargo (given by the 
proportion of out-sized cargo multiplied by the movement 
requirement for unit equipment). 


+ + X © © +*+ 


OUTREQMT(U,I,K) S$MOVEUE(U,I,K).. 
SUM ( (A,E1,E2,E3) $ ACCARGO(A, 'OUT'), 
SUM ( T $DSETX(U,A,I,K,El,E2,E3,T), TONSUE(U,A,1,K,E1l,E2,E3,T)) ) 
* UENOGO(U,I,K) sG2 CARGOP(U,'OUT') * MOVEUE(U,I,K) ; 


Total tonnage (in 100 tons) of unit equipment carried by over- 
sized cargo capable aircraft plus tonnage of unit equipment 

not delivered (UENOGO) must be greater than or equal to the 
movement requirement for over-sized and out-sized cargoes 
(given by the proportion of out-sized cargo + proportion of 
over-sized cargo multiplied by the movement requirement for 
unit equipment). This constraint is set up as such because an 
out-sized carrier can be used to carry over-sized cargo. 


+ + + X* f * ox 


OVERREQMT (U,I,K) SMOVEUE(U,I,K).. 
SUM ( (A,E1,E2,E3) $ACCARGO (A, *OVER'), 
SUM ( T $DSETX(U,A, I1,K,El,E2,E3,T), TONSUE(U,A,I,K,El,E2,E3., 110) 
+ UENOGO (U, I,K) 
=G= ( CARGOP(U,’OVER’) + CARGOP(U,’OUT’) ) * MOVEUE(U,I,K) ; 


* Total number (in 100s) of troops airlifted to the theater plus 
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* those not airlifted must be greater than or equal to the 
* movement requirement for troops. 


REQMTPAX(U,I,K) $MOVEPAX(U, I1,K).. 
SUM ( (A,E1,E2,E3), 
SUM ( T SDSETX(U, A.I K,E1,E2 E3,T), TPAX(U,A,1,K,El1,E2,E3,T) ) ) 
* PAXNOGO(U,I,K) =G= MOVEPAX(U,I,K) ; 


*++++++++++++++++ AIRCRAFT BALANCE CONSTRAINTS +++++++++++++++++ 


Aircraft balance constraint at origin airfields. On each day, 
the total number of aircraft assigned for airlift missions plus 
aircraft inventoried must be equal to the number of aircraft 
available from the last period plus new supply of a/c allocated 
to the base and aircraft that return from previous missions. 


+ + + + + 


ACBALI(A,I,T).. 
SUM ( (U,K,E1,E2,E3) $DSETX(U,A,1I,K,E1,E2,E3,T), 
X(U,A,1,K,E1,E2,E3,T) ) + H(A,1,T) 
=E= H(A,1,T-1)S(ORD(T) GT 1) + ALLOT(A,1I,T) SSUPPLY(A,T) + 
SUM ( (K,E1,E2,E3), 

SUM ( TP $( 
1SDSETY(A,1,K,E1,E2,E3,TP) AND 

15 ((ORD(TP) +ROUND (RCTTOI (A, I,K,E1,E2,E3)/24)) EQ ORD(T))), 
Y(A,I,K,E1,E2,E3,TP) ) ) ; 


* The sum of aircraft (newly made available) allocated to the 
different origin bases each day must be less than or equal 
* to the number of aircraft actually made available. 


* 


ACALLOT(A,T) SSUPPLY(A,T).. SUM(I, ALLOT(A,I,T)) =L= SUPPLY(A,T) ; 


Aircraft balance constraint at destination airfields. On each 
day, the number of each type of aircraft returning plus 

those that will return in the next period from a destination 
airfield must be equal to those waiting to return plus new 
arrivals at the destination. 


+ + + & * 


ACBALK(A,K,T).. 
SUNE (T El E2 E3) SDSETY (A, L K, El E2,E3;, T), 
Y(A,I,K,E1,E2,E3,T)) + HP(A,K,T) 
=E= HP(A,K,T-1)$(ORD(T) GT 1) + 
SUM ( (U,I,E1,E2,E3), 
SUM ( 
TP$( ((ORD(TP)+ROUND (CTTOK(A,I,K,E1,E2,E3)/24)) EQ ORD(T)) 
AND 1$DSETX(U,A,I,K,E1,E2,E3,TP)) , 
X(U,A,1,K,E1,E2,E3,TP) ) ) ; 


*++++++++++++ AIRCRAFT PHYSICAL LIMITATION CONSTRAINTS +++++++++++ 


Payload limitation: Total weight (equipment and men) carried 
by aircraft (A) on each route must be less than the critical 
payload for the route multiplied by the number of aircraft 
assigned for airlift missions. 


+ * + + 


ACWEIGHT(U,A,I,K,E1,E2,E3,T)SDSETX(U,A,I,K,E1,E2,E3,T) .. 
TONSUE (U,A, 1,K,E1,E2,E3,T) + PAXWT*TPAX (U,A, 1,K,E1,E2,E3,T) 
ZLE&MAXLOAD(A,I,K,EI,E2,E3) * X(U,A,I,K,El,E2,E3,T) ; 


* Cargo Space limitation: The total cargo space taken up by men 
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and equipment carried by aircraft (A) on each route must be 
less than the cargo space of the aircraft multiplied by the 
alrcraft loading efficiency and the number of aircraft 
assigned for airlift missions. 


+ +++ 


ACSPACE(U,A,I,K,E1,E2,E3,T) 
$(ACSQFT(A) $DSETX(U,A,I,K,E1,E2,E3,T) ).. 
( PAXSQFT(A) * TPAX(U,A,1,K,E1,E2,E3,T) ) + 
( UESQFT(U) * TONSUE(U,A,I,K,E1,E2,E3,T) ) 
=L= ACSQFT(A) * LOADEFF (A) * X(U,A,1,K,El,E2,E3,T) ; 


PAX carriage limitation: The total number of troops carried by 
aircraft (A) on each route must not be greater than the troop 
carriage capacity of the aircraft type multiplied by the number 
of aircraft assigned for airlift missions. 


+ + + + 


ACPAX(U,A,1,K,E1,E2,E3,T) S$DSETX(U,A,1,K,El, E2, B37 22 
TPAX(U,A,I,K,E1,E2,E3,T) =L= MAXPAX(A) * X(U,A,1I,K,E1,E2,E3,T) ; 


* Cargo space limitation. 


k+++++++++++++ AIRCRAFT UTILIZATION RATE CONSTRAINTS +++++++++++++ 


* The total number of flying hours (in 100s) consumed by each 

* aircraft type must be less than the utilization rate (in 100 hrs 

* per aircraft per day) times the number of aircraft initially 

* made available on day (T) times the remaining time periods where 
* these aircraft could be used. The productive time periods for 
* an aircraft made vailable on day (T) 1s equal to the total time 

* periods (CARD(T)) plus 1 minus ORD(T). For example, if the 

* total time periods of concern is 10, then an aircraft initially 
* made available on day 1 could be used for 10 days. 

ACURATE (A) .. 


SUM ( (U,I,K,E1,E2,E3,T) $DSETX(U,A,I,K,E1,E2,E3,T), 
X(U,A,1I,K,E1,E2,E3,T) * FLTTIME(A,I,K,E1,E2,E3) ) + 
SUM ( (I,K,E1,E2,E3,T) $DSETY(A,I,K,E1,E2,E3,T), 
Y(A,I,K,E1,E2,E3,T) * FLTTIME(A,I,K,E1,E2,E3) ) 
=L= SUM ( T $SUPPLY(A,T) , 
URATE (A) * ( CARD(T)+1-ORD(T) ) * SUPPLY(A,T) ); 


*++++ AIRCRAFT HANDLING CAPACITY OF AIRFIELD (MOG CONSTRAINT) ++++ 


These constraints model the throughput or handling capacity of 
an airfield; i.e., it will limit the number of aircraft using 
the airfield each day. As an airfield can serve as an origin, 
enroute, or destination, the total MOG consumed for each day 
would be a sum of individual consumptions. This total must not 
be greater than the MOG effciency factor multiplied by the MOG 
capacity (scaled to narrow body aircraft for comparison) of that 
airfield. A MOG efficiency factor is used as it is recognised 
that it is virtually impossible for an airfield to be running 
at full capacity for the whole day given the stochastic nature 
aircraft ground times. Unless aircraft can be sequenced to 
arrive at an airfield precisely one after another, congestion 
can occurs and the effective MOG utilization will be below the 
MOG capacity of the airfield. As an origin airfield, MOG 
consumed is equal to the total number of aircraft assigned for 
airlift missions scaled by the MOGREQ (normalize MOG 
requirement of aircraft) and the proportion of time (per day) 


* ox f + f f f f * f f + f ok ko 
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spent on a MOG spot.  Proportion of time spent is equal to the 
on-load time divided by 24 in this case. As an enroute 
airfield, the MOG consumed is equal to the total number of 
aircraft making enroute stops scaled by the MOGREQ factor and 
the proportion of time spent on a MOG spot (enroute ground time 
divided by 24). The MOG consumption when the airfield serves 
as a destination is computed in a similar manner except that 
the proportion of time spent is equal to the off-load time 
divided by 24. To prevent double counting of MOG consumption, 
conditional checks are done for enroute airfields. These are: 
1 Valid combination of indices (member of DSETX or DSETY) 
ii. Time of arrival is the period of concern 
iii. The aircraft came from somewhere (distance between the 
previous and current airfield is not equal to 0) 
iv. Time to enroute airfield is not equal to time to origin 
or destination airfield 


+ + + + + + +++ + ++ + +++ 


MOGUTILITY (AF,T).. 


SUM ( (U,A,K,E1,E2,E3) SDSETX(U,A,AF,K,E1,E2,E3,T) , 
X(U,A,AF,K,E1,E2,E3,T) * MOGREQ(AF,A) * GTIME(A,'ONLD') / 24 ) 
* MOG utilization when used as an origin airfield (I) 


+ 
SUM ( (U,A,I,K,E2,E3,TP) $ ( 
1SDSETX(U,A,I,K,AF,E2,E3,TP) AND 
15 (ORD(TP) +ROUND(CTTOE1 (A, I,K,AF,E2,E3)/24) EQ ORD(T)) AND 
1$DIST(I,AF) AND 
1S$ (CTTOE1(A,I,K,AF,E2,E3) NE CTTOK(A,I,K,AF,E2,E3)) ), 
X(U,A,I,K,AF,E2,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 


+ SUM ( (A,I,K,E2,E3,TP) S( 
1SDSETY (A,1I,K,AF,E2,E3,TP) AND 
1$ (ORD (TP) +ROUND (RCTTOE1 (A, 1,K,AF,E2,E3)/24) EQ ORD(T)) AND 
1SDIST(AF,E2) AND 
1S (RCTTOEl(A,I,K,AF,E2,E3) NE RCTTOI(A,I,K,AF,E2,E3)) ) , 
Y(A,I,K,AF,E2,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 
* MOG utilization when used as an El airfield 


+ 
SUM ( ({U,A,I,K,E1l1, E3, TP) ${ 
1SDSETX(U,A,I,K,E1,AF,E3,TP) AND 
1$ (ORD(TP) «*ROUND(CTTOE2(A,I,K,E1,AF,E3)/24) EQ ORD(T)) AND 
1SDIST(E1,AF) AND 
1$ (CTTOE2 (A,I,K,E1,AF,E3) NE CTTOK(A,1,K,E1,AF,E3)) ), 
X(U,A,I,K,ELl,AF,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 


* SUM ( (A,I,K,El,E3,TP) S( 
1SDSETY(A,I,K,El1,AF,E3,TP) AND 
1$ (ORD(TP) «*ROUND(RCTTOE2(A,I,K,E1,AF,E3)/24) EQ ORD(T)) AND 
1SDIST(AF,E3) AND 
1S(RCTTOE2 (A,I,K,E1,AF,E3) NE RCTTOI (A, 1,K,E1l,AF,E3)) ), 
Y(A,I,K,ELl,AF,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 
* MOG utilization when used as an E2 airfield 


+ 
SUM- (= (UA, I, K E1, E2 2) 91 
1$DSETX (U,A,1,K,E1,E2,AF,TP) AND 
1$ (ORD (TP) +ROUND (CTTOE3 (A, 1I,K,E1,E2,AF) /24) EQ ORD(T)) AND 
1SDIST(E2,AF) AND 
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SUM ( 


1SDIST(AF,K) AND 


MODEL AIRLIFT  /ALL/; 


Include File Summary 


GLOBAL TYPE LOCAL 
0 INPUT 0 
24 INCLUDE 24 
49 INCLUDE Or 
61 INCLUDE 30 
69 INCLUDE 33 
89 INCLUDE 47 
107 INCLUDE 59 
130 INCLUDE 62 
139 INCLUDE 65 
147 INCLUDE 68 
162 INCLUDE 71 
175 INCLUDE 74 
190 INCLUDE 81 
237 INCLUDE 84 
267 INCLUDE 91 
298 INCLUDE 98 
323 INCLUDE 102 
350 INCLUDE 105 
375 INCLUDE 108 
400 INCLUDE 111 
407 INCLUDE 115 
439 INCLUDE 123 
445 INCLUDE 128 
478 INCLUDE 138 
507 INCLUDE 142 
520 INCLUDE 146 
533 INCLUDE 150 
547 INCLUDE 154 
564 INCLUDE 158 
579 INCLUDE 162 


1$ (CTTOE3 (A, 1,K,E1,E2,AF) NE CTTOK(A,1,K,E1,E2,AF)) Kip 
X(U,A,1I,K,E1,E2,AF,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 


(A,I,K,E1,E2,TP) $( 
1$DSETY (A, 1,K,E1,E2,AF,TP) AND 
1$ (ORD(TP) +ROUND(RCTTOE3 (A, 1,K,E1,E2,AF)/24) EQ ORD(T)) AND 


1$ (RCTTOE3 (A,I,K,E1,E2,AF) NE RCTTOI(A,1I,K,E1,E2,AF)) ), 
Y(A,I,K,E1,E2,AF,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 
MOG utilization when used as an E3 airfield 


(U,A,I,E1,E2,E3,TP) $( 

1$DSETX(U,A,I,AF,E1,E2,E3,TP) AND 

1$ (ORD(TP) +ROUND(CTTOK(A,1,AF,E1,E2,E3)/24) EQ ORD(T)) ), 
X(U,A,1I,AF,E1,E2,E3,TP) * MOGREQ(AF,A) * GTIME(A,’OFFLD’) / 24) 
MOG utilization when used as a destination airfield (K) 


=L= MOGEFF * MOGCAP(AF,'NB') ; 


SOLVE AIRLIFT using RMIP minimizing Z ; 


FILE NAME 


/home/limt/airlift6.gms 
./home/limt/unitname.dat 
./home/limt/acname.dat 

. /home/limt/craf.dat 
./home/limt/milac.dat 
./home/limt/periods.dat 
./home/limt/afname.dat 
./home/1limt/orgname.dat 
./home/limt/destname.dat 
./home/1limt/einame.dat 
./home/limt/e2name.dat 
./home/limt/e3name.dat 
./home/limt/vroutex.set 
./home/limt/vroutey.set 
./home/limt/moveue.dat 

. /home/limt/cargop.dat 
./home/limt/movepax. dat 
./home/limt/ald.dat 
./home/limt/rdd.dat 
./home/limt/maxlate.dat 
./home/limt/latepen.dat 
./home/1limt/nogopen.dat 
./home/limt/uesqft.dat 

. /home/limt/acsupply.dat 
./home/limt/acsize.dat 
./home/limt/accargo.dat 
. /home/limt/acmaxpax.dat 
. /home/limt/paxsqft.dat 
./home/limt/acsqft.dat 
./home/limt/loadeff.dat 
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595 INCLUDE 
608 INCLUDE 
621 INCLUDE 
644 INCLUDE 
656 INCLUDE 
682 INCLUDE 
709 INCLUDE 
735 INCLUDE 


COMPILATION TIME 
Model Statistics 
MODEL STATISTICS 
BLOCKS OF EQUATIONS 
BLOCKS OF VARIABLES 


NON ZERO ELEMENTS 


GENERATION TIME 


EXECUTION TIME 


STEP SUMMARY: 0. 
.680 COMPILATION 
.700 EXECUTION 
.090 CLOSEDOWN 
.480 TOTAL SECONDS 


44 


45 


Solution Report 


166 ./home/limt/acspeed.dat 
170 ./home/limt/gtime.dat 
175 ./home/limt/urate.dat 
188 ./home/limt/acrange.dat 
192 ./home/limt/acload.dat 
209 ./home/limt/afcoord.dat 
218 ./home/limt/mogcap.dat 
226 ./home/limt/mogeff.dat 


^ 0.680 SECONDS VERID AIX-00-064 


SOLVE AIRLIFT USING RMIP FROM LINE 1218 


1:3 SINGLE EQUATIONS 6349 
10 SINGLE VARIABLES 8723 
38614 DISCRETE VARIABLES 5500 


42.590 SECONDS 


- 44.700 SECONDS VERID AIX-00-064 


010 STARTUP 


SOLVE AIRLIFT USING RMIP FROM LINE 1218 


SOLVE SUMMARY 


MODEL AIRLIFT 


TYPE RMIP 
SOLVER OSL 


τν. SOLVER STATUS 


OBJECTIVE 2 
DIRECTION MINIMIZE 
FROM LINE 1218 


1 NORMAL COMPLETION 


**** MODEL STATUS 1 OPTIMAL 

**x** OBJECTIVE VALUE 37.0139 

RESOURCE USAGE, LIMIT 28.690 18000.000 

ITERATION COUNT, LIMIT 3001 400000 

OSL Release 2, GAMS Link level 3 --- AIX RS/6000 1.3.045-017 
Work space allocated -- 5.54 Mb 

**** REPORT SUMMARY 0 ΝΟΝΟΡΤ 


EXECUTION TIME 


0 INFEASIBLE 
0 UNBOUNDED 


= 0.230 SECONDS VERID AIX-00-064 
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USER: Operations Research Department 
Naval Postgraduate School 


«eee FILE SUMMARY 


INPUT /home/limt/airlift6.gms 
OUTPUT /home/limt/airlift6.1lst 
SAVE /home/limt/final.go? 


**** WARNING - COMPILER OPTIONS ARE NON DEFAULT 
INLINECOM { 

STEP SUMMARY: .180 STARTUP 

.000 COMPILATION 

.230 EXECUTION 

.150 CLOSEDOWN 

.560 TOTAL SECONDS 


ΟΟΟΟΟ 
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G940405-1442AX- AIX 


1215 
1220 
1221 
1222 
1223 
1224 


1225 
1226 
1227 
1228 
1229 
1230 
1231 
1232 


1233 
1234 
1235 
1236 
1237 
1238 
1239 
1240 
1241 


1242 
1243 
1244 
1245 
1246 
1247 
1248 
1249 


1250 
1251 
1252 
1253 
1254 
1255 
1256 
1257 
1258 
1259 
1260 
1261 
1262 
1263 
1264 
1265 
1266 
1267 
1268 
1269 
1270 
1271 
1272 
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*tt4tt 4 xGBRG AÁG MG AG MX M — 44484444 SUMMARY REPORTS «*tt644ttt4646 6644 t6 t. 444 t 44 
f+++++++++++++++++++++ DEMAND SATISFACTION +++++++++++++++++++++++ 


PARAMETER TONSONT(U,I,K) unit U's equipment (in 100 tons) 
* delivered on time; 
TONSONT(U,I,K) SMOVEUE(U,I,K) 
Ba SUM( (A,EL,E2,E3,T), 
TONSUE.L(U,A,I,K,E1,E2,E3,T) 
$((ORD(T) *ROUND(CTTOK(A,I,K,E1,E2,E3)/24) LE RDD(U)) 
SDSETX(U,A,I,K,E1,E2,E3,T)) 
js 


PARAMETER TONSLATE (U, I, K) unit U's equipment (in 100 tons) 
la delivered late; 
TONSLATE (U,I,K) SMOVEUE(U,I,K) 
= SUM( (A,E1,E2,E3,T), 
TONSUE.L(U,A,I,K,E1,E2,E3,T) 
S ((ORD(T) +ROUND(CTTOK(A,1,K,E1,E2,E3)/24) GT RDD(U)) 
SDSETX(U,A,I,K,E1,E2,E3,T)) 
)5 


PARAMETER PAXONT(U,I,K) number (in 100s) of unit U's troops 
* airlifted on time; 
PAXONT(U,I,K) SMOVEPAX(U,I,K) 
= SUM( (A,E1,E2,E3,T), 
TPAX.L(U,A,1I,K,E1,E2,E3,T) 
S$ ( (ORD (T) +ROUND (CTTOK(A,1I,K,E1,E2,E3) /24) LE RDD(U) ) 
SDSETX(U,A,I,K,E1,E2,E3,T)) 
E 


PARAMETER PAXLATE (U,I,K) number (in 100s) of unit U's troops 
* airlifted late; 
PAXLATE(U,I,K) SMOVEPAX(U,I,K) 
= SUM( (A,E1,E2,E3,T), 
TPAX.L(U,A,I,K,E1,E2,E3,T) 
$ ( (ORD (T) +ROUND (CTTOK(A,1I,K,E1,E2,E3) /24) GT RDD(U) ) 
SDSETX(U,A,I,K,E1,E2,E3,T)) 
J5 


OPTION TONSONT:2:0:3; 
DISPLAY TONSONT; 
OPTION TONSLATE:2:0:3; 
DISPLAY TONSLATE; 
OPTION UENOGO:2:0:3; 
DISPLAY UENOGO.L; 
OPTION PAXONT:2:0:3; 
DISPLAY PAXONT; 
OPTION PAXLATE:2:0:3; 
DISPLAY PAXLATE; 
OPTION PAXNOGO:2:0:3; 
DISPLAY PAXNOGO.L; 


k++++++++++++++++++ AIRCRAFT UTILIZATION RATE ++++++++++++++++++++ 
PARAMETER ACTUALUR(A) actual utilization rate (in hrs per day per 


* aircraft); 
ACTUALUR(A) = 
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( SUM ( (U,I,K,E1,E2,E3,T) S$SDSETX(U,A,1,K,E1,F2,E3 DP 
X.L(U,A,I1,K,EB1,E2,E3,T) * FLTTIME(A,I,K,EËE1,F27/E3 PK 
SUM ( (I,K,E1,E2,E3,T) $DSETY(A,I1,K,E1,E2,E3,"T) , 
Y.L(A,1,K,E1,E2,E3,T) * FLTTIME(A,1,K,E1,E2,E3) ) ) * 100 
/ SUM ( T, ( CARD(T)+1-ORD(T) ) * SUPPLY(A,T) SSUPPLY(A,T) ); 
DISPLAY ACTUALUR; 


Y ++T+T++++++++++++++++++++ MOG CONSUMPTION +++++++++++++++++++++++++ 


Reports the MOG consumption of an airfield. If you set the MOG 
efficiency factor to 0.9, and the percentage of MOG capacity 
used, MOGUSED(AF,T), is 90, then it means that the airfield has 
"MOGed" out; i.e., the airfield's aircraft handling capacity is 
a binding constraint. 


+t f X co 


PARAMETER MOGUSED (AF, T) percentage of MOG capacity used ; 
MOGUSED(AF,T) = ( 
SUM ( (U,A,K,El,E2,E3) $DSETX(U,A,AF,K,E1,E2,E3,T) , 
X.L(U,A,AF,K,E1,E2,E3,T) * MOGREQ(AF,A) * GTIME(A,’ONLD’) / 24 ) 
* MOG utilization when used as an origin airfield (I) 


+ 

SUM ( (U,A,I,K,E2,E3,TP) $( 
1$DSETX (U,A, 1,K,AF,E2,E3,TP) AND 
1$ (ORD (TP) +ROUND(CTTOE1 (A, 1,K,AF,E2,E3)/24) EQ ORD(T)) AND 
1$DIST(I,AF) AND 
1$ (CTTOE1(A,1,K,AF,E2,E3) NE CTTOK(A, 1,K,AF,E2, ES) eee 
X.L(U,A,I,K,AF,E2,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 


+ SUM ( (A,1I,K,E2,E3,TP) $( 
1SDSETY(A,I,K,AF,E2,E3,TP) AND 
1$ (ORD (TP) +ROUND (RCTTOE1(A,I,K,AF,E2,E3) /24) EQ ORD(T)) AND 
1$DIST(AF,E2) AND 
1S(RCTTOE1(A,1,K,AF,E2,E3) NE RCTTOI(A,I,K,AF,E2,E3)) ), 
Y.L(A,I,K,AF,E2,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 
* MOG utilization when used as an El airfield 


+ 

SUM ( (U,A,I,K,E1,E3,TP) $( 
1$DSETX(U,A,I,K,E1,AF,E3,TP) AND 
1$ (ORD (TP) +ROUND (CTTOE2 (A, 1,K,E1,AF,E3)/24) EQ ORD(T)) AND 
1$DIST(E1,AF) AND 
1$ (CTTOE2(A,I,K,El1,AF,E3) NE CTTOK(A,I,K,E1,AF,E3)) ), 
X.L(U,A,1I,K,E1,AF,E3,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 


+ SUM ( (A,I,K,El1,E3,TP) S( 
l1SDSETY(A,I,K,El,AF,E3,TP) AND 
1$ (ORD(TP) «ROUND (RCTTOE2(A,I,K,E1,AF,E3)/24) EQ ORD(T)) AND 
1$DIST(AF,E3) AND 
1$ (RCTTOE2 (A,I,K,E1,AF,E3) NE RCTTOI(A,I,K,E1,AF,E3)) ), 
Y.L(A,I,K,E1,AF,E3,TP) * MOGREQ(AF,A) * GTIME(A,’ENR’) / 24 ) 
* MOG utilization when used as an E2 airfield 


+ 

SUM ( (U,A,I,K,E1,E2,TP) S( 
1$DSETX(U,A,I,K,E1,E2,AF,TP) AND 
1$ (ORD (TP) +ROUND(CTTOE3 (A, 1,K,E1,E2,AF)/24) EQ ORD(T)) AND 
1$DIST(E2,AF) AND 
1$ (CTTOE3 (A,I,K,E1,E2,AF) NE CTTOK(A,I,K,E1,E2,AF)) ), 
X.L(U,A,I,K,E1,E2,AF,TP) * MOGREQ(AF,A) * GTIME(A,’ENR’) / 24 ) 
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f 


1333 


1334 + SUM ( (A,I,K,E1,E2,TP) $( 

1335 l$DSETY(A,I,K,E1,E2,AF,TP) AND 

1336 1$ (ORD (TP) +ROUND(RCTTOE3 (A, 1,K,E1,E2,AF)/24) EQ ORD(T)) AND 

1337 1SDIST(AF,K) AND 

1338 1$ (RCTTOE3 (A,I,K,E1,E2,AF) NE RCTTOI(A,1I,K,E1,E2,AF)) ), 

1339 Y.L(A,I,K,ELl,E2,AF,TP) * MOGREQ(AF,A) * GTIME(A,'ENR') / 24 ) 

1340 * MOG utilization when used as an E3 airfield 

1341 

1342 + 

1343 SUM ( (U,A,I,El,E2,E3,TP) $( 

1344 l1$DSETX(U,A,I,AF,E1,E2,E3,TP) AND 

1345 1$ (ORD (TP) +ROUND (CTTOK(A,I,AF,E1,E2,E3)/24) EQ ORD(T)) ), 

1346 X.L(U,A,I,AF,E1,E2,E3,TP) *MOGREQ(AF,A) *GTIME(A,'OFFLD') / 24) ) 

1347 * MOG utilization when used as a destination airfield (K) 

1348 

1349 / MOGCAP(AF,'NB') * 100 ; 

1350 

1351 OPTION MOGUSED:2:0:4; 

1352 DISPLAY MOGUSED; 

1353 

1354 *++++++++++ AIRCRAFT ALLOCATION AND MISSION ASSIGNMENTS) +4+++4+4++++++ 

1355 

1356 OPTION ALLOT:2:0:2; 

1357 DISPLAY ALLOT.L; 

1358 OPTION X:2:0:1; 

1359 DISPLAY X.L; 

1360 OPTION Y:2:0:1; 

1361 DISPLAY Y.L; 

1362 OPTION TONSUE:2:0:1; 

1363 DISPLAY TONSUE.L; 

1364 OPTION TPAX:2:0:1; 

1365 DISPLAY TPAX.L; 

1366 

COMPILATION TIME Ξ 0.060 SECONDS VERID AIX-00-064 

---- 1258 PARAMETER TONSONT unit U's equipment (in 100 tons) 
delivered on time 

UNITA.XDAT.FFTJ 79.02, UNITB.TYFR.FFTJ 6.92, UNITD.TMKH.FFTJ 114.94 

UNITE.TYFR.UGZX 27.96, UNITF.TMKH.UGZX 183.96, UNITG.TYFR.UGZX 5.03 

UNITH.XDAT.UGZX 3.01, UNITL.TYFR.UGZX 4.91, UNITM.TYFR.UGZX 3.24 

UNITN.TYFR.UGZX 5.03, UNITO.TYFR.PKVV 3.07, UNITP.TYFR.PKVV 5 2c 

UNITQ.TMKH.PKVV 959. UNITR.TMKH.UGZX 10.13, UNITS.TMKH.UGZX 10.13 

UNITT.XDAT.UGZX 4.78 

---- 1260 PARAMETER TONSLATE unit U's equipment (in 100 tons) 
delivered late 

UNITC.XDAT.FFTJ 42.35, UNITF.TMKH.UGZX 54.31, UNITH.XDAT.UGZX 7.12 

UNITI.QFQE.UGZX 5.03, UNITJ.XDAT.UGZX 10.13, UNITK.TYFR.UGZX 5.03 

UNITM.TYFR.UGZX 1.98, UNITO.TYFR.PKVV 1.84, UNITQ.TMKH.PKVV 0.54 

UNITT.XDAT.UGZX 5.35 

---- 1262 VARIABLE UENOGO.L unit U's equipment (in 100 tons) that 
are not airlifted 

UNITA.XDAT.FFTJ 107.75, UNITC.XDAT.FFTJ 20.88, UNITE.TYFR.UGZX 635.73 
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UNITA. 
UNITD. 
UNITG. 
UNITJ. 
UNITM. 
UNITP. 


UNITS 


UNITC. 


ES 
747C 


QFQE. 
QFOE. 
QFQE. 
QFQE. 
OEOE: 
OQFOE. 
QFQE. 
TMKH. 
TMKH. 
TMKH. 
ΤΜΚΗ. 
TMKH. 
TYFR. 
TYFR. 
TIER 
TYER. 
TYFR. 
TIPR. 
XDAT. 
XDAT 
XDAT. 
XDAT. 
XDAT. 
PIE 
FFTJ. 
FETI. 
τα. 
FETI. 
UGZX. 
UGZX. 


1264 PARAMETER PAXONT 


XDAT 


1266 PARAMETER PAXLATE 


XDAT. 
TMKH. 
TIER? 


BIER. 
TIERS 
. TMKH. 


FFTJ 126.98, 
FFTJ 97.76, 
UGZX 379725 
. UGZX 4.35, 
UGZX 3972. 
PKVV 3.721 
UGZX 8.46, 


XDAT.FFTJ 22.16, 


UNITB.TYER.FEFTJ "13. 
UNITE.TYFR.UGZX 199. 


number (in 100s) of unit U's troops 
airlifted on time 


UNITH.XDAT.UGZX 
UNITK.TYFR.UGZX 
UNITN.TYFR.UGZX 
UNITQ.TMKH.PKVV 
UNITT.XDAT.UGZX 


oo OO Ú) UJ OD 


UNITC.XDAT.FFTJ 30.43 


UNITF.TMKH.UGZX 190. 


UNITI.QFQE.UGZX 
UNITL.TYFR.UGZX 
UNITO.TYFR.PKVV 
UNITR.TMKH.UGZX 


Oh bh W 
~J 
> 


number (in 100s) of unit U’s troops 
airlifted late 


UNITJ.XDAT.UGZX 4.11 


1268 VARIABLE  PAXNOGO.L 


( ALL 


1279 PARAMETER ACTUALUR 


9.32, 
3.25 


1352 PARAMETER MOGUSED 


DAY2 
DAY6 
DAY11 
DAY15 
DAY19 
DAY23 
DAY27 
DAY5 
DAY12 
DAY16 
DAY20 
DAY24 
DAY1 
DAY7 
DAY12 
DAY16 
DAY20 
DAY24 
DAY1 


.DAY7 


DAY11 
DAY15 
DAY20 
DAY11 
DAY15 
DAY19 
DAY23 
DAY27 
DAY3 

DAY7 


C17 


9 39, 


DC10 0.44 


19.48, 
70:60, 
90:007 
ο οὐ 
20.00; 
90 
90 
ο ο 
13.83. 
34.18, 
30995 
49 
2 55, 
24.30, 
82.42, 
76 
16 
28 
5.04, 
10.40, 
42 


.64, 


31 
0.46, 
16 
4.93, 
75 
85. 
8.88, 
90. 


ο 
ο 


X15 


.54, 
. 105 
720, 


19:07, 
2307 


QFQE 
OFOE. 
QFQE 
QFQE 
QFQE 
QFQE. 
QFQE. 
TMKH 
TMKH. 
TMKH. 
TMKH. 
TMKH. 
TIER. 
TYFR. 
TYFR. 
TER: 
TILER: 
TYFR. 
XDAT. 
XDAT. 
XDAT. 
XDAT. 
EFIJ. 
FFTJ. 
FFTJ. 
FFTJ. 
FPIJ. 
PRI: 
UGZX. 
UGZX 


C141 3.05, 


.DAY3 


DAY7 


.DAY12 
.DAY16 
.DAY20 


DAY24 
DAY28 


.DAY6 


DAY13 
DAY17 
DAY21 
DAY25 
DAY2 

DAY8 

DAY13 
DAY17 
DAY21 
DAY25 
DAY2 

DAY8 

DAY12 
DAY16 
DAY7 

DAY12 
DAY16 
DAY20 
DAY24 
DAY28 
DAY4 


.DAY8 


number (in 100s) of unit U's troops 
not airlifted 


0.00 


) 


actual utilization rate 
day per aircraft) 


C130 3.793, 


747P 


1.367 


percentage of MOG capacity 
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QFOE. 
Ες. 
SCENES 
QFQE. 
QFQE. 
QFOE. 
OFOE- 
TMKH 
TMKH. 
TMKH. 
TMKH. 
TMKH. 
TIER 
τινα. 
DERE 
τις. 
TYER. 
TYE Re 
XDAT . 
XDAT. 
XDAT 
XDAT 
IPOD 
ELE 
PELJ. 
FFTJ. 
FFTJ. 
UGZX. 
UGZX. 
UGZX. 


DAY4 

DAY8 

DAY13 
DAY17 
DAY21 
DAY25 
DAY29 


. DAY7 


DAY14 
DAY18 
DAY22 
DAY26 
DAYS 

DAY10 
DAY14 
DAY18 
DAY22 
DAY26 
DAY3 

DAY9 


.DAY13 
.DAY17 


DAY8 
DAY13 
DAY17 
DAY21 
DAY25 
DAY1 
DAY5 
DAY9 


.40, 
. O87 
οο 
200 
ος, 
EO 
.84, 
aT Sy 
27097 
207; 


QFQE 
QFQE 


QFQE. 
.DAY18 
.DAY22 
.DAY26 
TMKH. 
TMKH. 
TMKH. 
TMKH. 
.DAY23 
TMKH. 
TYFR 
TYFR: 
TYERS 
TYFR. 
.DAY23 
TYER 
XDAT. 
XDAT. 
XDAT. 
XDAT. 
.DAY9 
FFTJ. 
FFTJ. 
FFTJ. 
FFIJ. 
.DAY2 
UGZX. 
.DAY10 


QFQE 
QFQE 
QFQE 


TMKH 


IE 


FFTJ 


UGZX 


UGZX 


.DAY5 
. DAY9 


DAY14 


DAY4 

DAY11 
DAY15 
DAY19 


DAY27 
DAY6 

DAY11 
DAY15 
DAY19 


DAY27 
DAY6 

DAY10 
DAY14 
DAY19 


DAY14 
DAY18 
DAY22 
DAY26 


DAY6 


(in hrs per 


.DAY14 
.DAY18 
.DAY22 
. DAY6 

.DAY18 
. DAYS 

.DAY15 
.DAY20 
.DAY25 
.DAY9 

. DAY14 
.DAY20 
.DAY10 
.DAY16 
. DAY7 

DAY11 
.DAY16 
.DAY20 
. DAY9 

.DAY21 


DAY9 


. DAY14 


DAY18 


UGZX.DAY11 90.00, UGZX.DAY12 90.00, UGZX.DAY13 90.00, UGZX 
UGZX.DAY15 90.00, UGZX.DAY16 90.00, UGZX.DAY17 90.00, UGZX 
UGZX.DAY19 90.00, UGZX.DAY20 90.00, UGZX.DAY21 90.00, UGZX 
UGZX.DAY26 5.23, UGZX.DAY28 31.35, PKVV.DAYS 90.00, PKVV 
PKVV.DAY7 90.00, PKVV.DAY8 3.28,  KNMD.DAY7 0.73, KNMD 
KNMD.DAY19 6.00, FXSB.DAY7 4.03, FXSB.DAY8 7.50, FXSB 
FXSB.DAY11 4.03, WWYK.DAY13 4.69, WWYK.DAY14 43.74, WWYK 
WWYK.DAY17 53.02, WWYK.DAY18 50.16, WWYK.DAY19 50.01, WWYK 
WWYK.DAY21 90.00, WWYK.DAY22 9.37, WWYK.DAY24 56.25, WWYK 
NNGX.DAY6 11.45, NNGX.DAY7 12.36, NNGX.DAY8 11.45, NNGX 
NNGX.DAY10 1.35, NNGX.DAY11 25.63, NNGX.DAY13 25.63, NNGX 
NNGX.DAY15 18.75, NNGX.DAY16 59.42, NNGX.DAY17 43.68, NNGX 
NNGX.DAY22 51.43, UTKY.DAY?7 4.61, UTKY.DAYS8 2.88, UTKY 
ΠΠ ΡΛΥΙ3 25.71,  UTKY.DAY14 21.11, UTKY.DAY15 25.71, UTKY 
UTKY.DAY17 25.71, UTKT.DAY20;25.71, .UTKY.DAY21 12.86, ZNRE 
ZNRE.DAY8 8.02, ZNRE.DAYS 6.45, ZNRE.DAY10 0.94, ZNRE 
ZNRE.DAY13 21.22, ZNRE.DAY14 14.78, ZNRE.DAY15 18.00, ZNRE 
ZNRE.DAY17 27.00, ZNRE.DAY18 9.00, ZNRE.DAY19 18.00, ZNRE 
ZNRE.DAY21 27.00, HTDS.DAY7 20.01, HTDS.DAY8 10.75, HTDS 
HTDS.DAY11 3.14, HTDS.DAY12 59.80, HTDS.DAY18 90.00, HTDS 
FGDC.DAY6 43.25, FGDC.DAY7 40.08, FGDC.DAY8 3.24, FGDC. 
FGDC.DAY10 16.12, FGDC.DAY12 90.00, FGDC.DAY13 90.00, FGDC 
FGDC.DAY15 90.00, FGDC.DAY16 90.00, FGDC.DAY17 90.00, FGDC. 
FGDC.DAY19 90.00, FGDC.DAY20 90.00, FGDC.DAY21 90.00 
---- 1357 VARIABLE ALLOT.L new aircraft assets made available 
allotted to base (I). 
C5 . TYFR.DAY1 1.44, C5 .XDAT.DAY1 8.56 
Gs .XDAT.DAY6 30.00, ος .XDAT.DAY11 40.00 
C17 .TYFR.DAY1 2 00; C17 .TYFR.DAY6 6.00 
C17 .TYFR.DAY11 8.00, C141.QFQE.DAY1 3.29 
C141.QFQE.DAY6 5 C141.TYFR.DAY1 5.02 
C141.TYFR.DAY6 53 B3. C141.XDAT.DAY1 11.69 
C141.XDAT.DAY11 80.00, C130.TYFR.DAY1 20.00 
C130.TYFR.DAY6 60.00, C130.TYFR.DAY11 80.00 
747P.TYFR.DAY6 5 51, 747P.TYFR.DAY11 30.00 
747P.TYFR.DAY16 20.00, 747P.XDAT.DAY6 3.49 
747C.QFQE.DAY6 2.68, 747C.TYFR.DAY6 1.19 
747C.TYFR.DAY11 7.34, 747C.TYFR.DAY16 8.37 
747C.XDAT.DAY6 1.1.3: 747C.XDAT.DAY11 17.66 
747C.XDAT.DAY16 29.35, DC10.TMKH.DAY11 JO 
DC10 .TYFR.DAY11 17.90, DC10.TYFR.DAY16 6.29 
DC10.XDAT.DAY6 5.00; DC10.XDAT.DAY16 9.67 
---- 1359 VARIABLE X.L number of airlift missions (aircraft) 
UNITA.C141.XDAT.FFTJ.QFQE.FFTJ.FFTJ.DAY12 11-55 
UNITA.C141.XDAT.FFTJ.QFQE.FFTJ.FFTJ.DAY13 9.98 
UNITA.C141.XDAT.FFTJ.QFQE.FFTJ.FFTJ.DAY14 8.48 
UNITA.C141.XDAT.FFTJ.QFQE.FFTJ.FFTJ.DAY15 33.48 
UNITA.C141.XDAT.FFTJ.QFQE.FFTJ.FFTJ.DAY16 3.51 
UNITA.747P.XDAT.FFTJ.XDAT.ZNRE.FGDC.DAY12 3.87 
UNITA.747P.XDAT.FFTJ.XDAT.ZNRE.FGDC.DAY16 10.80 
UNITA.747P.XDAT.FFTJ.NNGX.TYFR.FFTJ.DAY16 235.98 
UNITA.747C.XDAT.FFTJ.XDAT.ZNRE.FGDC.DAY19 10.80 
UNITA.747C.XDAT.FFTJ.XDAT.ZNRE.FGDC.DAY20 21.60 
UNITA.747C.XDAT.FFTJ.NNGX.TYFR.FFTJ.DAY13 21.53 
UNITA.747C.XDAT.FFTJ.NNGX.TYFR.FFTJ.DAY14 8.47 
UNITA.747C.XDAT.FFTJ.NNGX.TYFR.FFTJ.DAY15 15.75 


19 


QFQE. 


.DAY16 
.DAY17 
.DAY17 
.DAY20 
.DAY17 
.DAY16 
.DAY18 
.DAY17 
.DAY19 
.DAY11 
.DAY11 
.DAY6 

. DAY? 

.DAY8 

.DAY10 
.DAY11 
.DAY6 

.DAY8 

.DAY21 
.DAY22 
.DAY24 
.DAY25 
.DAY27 
.DAY22 
.DAY26 
.DAY22 
.DAY23 
. DAY24 
.DAY25 
.DAY26 
.DAY22 
.DAY27 
.DAY21 
. DAY22 
.DAY26 
.DAY11 
.DAY12 
.DAY11 
.DAY13 
.DAY16 
.DAY11 
.DAY16 
.DAY11 
.DAY11 
.DAY12 
.DAY13 
.DAY14 
.DAY15 
.DAY16 
.DAY17 
.DAY18 
.DAY19 
.DAY20 
.DAY21 
.DAY11 
.DAY15 
: .DAY17 
. UGZX. 
. UGZX. 
. UGZX. 


DAY19 
DAY20 
DAY11 


80 


UNITT. 


nol 7 


ΠΕΠ. 


εκ. 
ΕΞ. 
ER. 
T PR. 
ΓΤ. 
ΕΕ. 


ΠΙΕΡ. 


JQFQE. 
.QFQE 
.QFQE 
.QFQE 
.QFQE 
.QFQE 
.QFQE 
. QFQE 
. QFQE 
OIYER 
oo ER. 
.QFQE 
Ρος 
ΠΕΠ. 
.QFQE 
.NNGX. 
.NNGX. 


UGZX. 
.UGZX. 
.UGZX. 
. UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
τες. 
απ κ 
ΠΕΡ. 


STIER: 
Εκ. 
ΓΕΝ. 
.UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
. UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
“TIER: 
TIPR. 
το κ. 
HTDS. 
ot VER 
οκ: 
SUI REOS 
EYER: 
.TYFR. 
“τος. 
τες. 
ΕΕ. 
ΕΕ. 
τες, 
ot GER 
ΕΙ. 
TAER: 
JEVER: 
.IYER. 
UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
.UGZX. 
. UGZX. 
.UGZX. 
. UGZX. 
. UGZX. 
. UGZX. 
UGZX. 
.UGZX. 
.UGZX. 
UGZX. 
.UGZX. 


TYFR. 


4DAYI2 
.DAY13 
.DAY14 
.DAY16 
.DAY17 
.DAY18 
.DAY21 
.DAY20 
.DAY12 
.DAY12 
.DAY16 
. DAY1 
. DAY1 
. DAY 2 
. DAY1 
. DAY1 
.DAY6 
. DAY 2 
. DAY 3 
. DAY2 
.DAY6 
.DAY6 
.DAY6 
.DAY2 
.DAY6 
.DAY6 
. DAY7 
. DAY8 
.DAY2 
.DAY6 
.DAY6 
.DAY6 
.DAY6 
. DAY6 
. DAY6 
. DAY6 
.DAY6 
.DAY6 
.DAY6 
.DAY6 
.DAY6 
.DAY6 
.DAY6 
.DAY4 
.DAY4 
.DAY6 
.DAY7 
.DAY4 
.DAY5 
.DAY6 
.DAY4 
.DAY4 
.DAY6 
.DAY6 


DAY7 


. DAY7 
. DAY7 
. DAY7 
.DAY7 
. DAY7 


81 


μι P 


N 


pp = 
PP P O Q BP Ñ BP BP jJ O Ñ O Q Ñ P. i£ O O P Ñ P @ Q Q) +J O O O @ G mo mo 


N HP 
Ποσο” ο ο «Οδ ο”, ΟΝ 


= 


ω 


w P BP mpm 
Pb PÑ OP OPF ο)” 


87 
76 


UNITT.747C.XDAT.UGZX.NNGX.TYFR.UGZX.DAY9 3 
UNITT.747C.XDAT.UGZX.NNGX.HTDS.UGZX.DAY10 ο 


1361 VARIABLE Y.L number of type A aircraft recovered 


CS .TMKH.FFTJ.QFOE.FETJ EET 'DAY26 0 59 
C5 .TMKH.FFTJ.TYFR.FFTJ.FFTJ.DAY12 34.64 
CS .TMKH.FFTJ .WWYK.FFTJ.FFTJ.DAY18 5.36 
CS .TMKH.FFTJ.WWYK.FFTJ.FFTJ.DAY23 48.00 
C5 .TMKH.UGZX.QFQE.UGZX.UGZX.DAY3 10.00 
C5 .TMKH.UGZX.QFQE.UGZX.UGZX.DAY9 37.92 
C5  .TMKH.UGZX.QFQE.UGZX.UGZX.DAY15 24.90 
C5  .TMKH.UGZX.QFQE.UGZX.UGZX.DAY21 25.44 
CS  .TMKH.UGZX.QFQE.UGZX.UGZX.DAY23 9.22 
C5  .TMKH.UGZX.TYFR.UGZX.UGZX.DAYS 3 ον 
C5  .TMKH.UGZX.TYFR.UGZX.UGZX.DAY6 4.89 
CS  .TMKH.UGZX.TYFR.UGZX.UGZX.DAY12 2.68 
C5 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY15 4.50 
C5 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY20 5.15 
C5 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY23 22.78 
C5  .TMKH.UGZX.WWYK.UGZX.UGZX.DAY13 37.32 
CS  .TMKH.UGZX.WWYK.UGZX.UGZX.DAY14 13.94 
CS  .TMKH.UGZX.WWYK.UGZX.UGZX.DAY16 45.24 
C5  .TMKH.UGZX.WWYK.UGZX.UGZX.DAY17 29.40 
C5  .TMKH.UGZX.WWYK.UGZX.UGZX.DAY18 37.32 
ES TMKH.UGZX.WWYK.UGZX.UGZX.DAY19 20.08 
ES TMKH.UGZX.WWYK.UGZX.UGZX.DAY20 49.41 
C5 .TMKH.PKVV.TYFR.UGZX.PKVV.DAYS 2.08 
C17 .TMKH.FFTIJ.WWYK.FFTJ.FFTJ.DAY24 167098 
C17 .TMKH.UGZX.QFQE.UGZX.UGZX.DAY3 2.00 
C17 .TMKH.UGZX.QFQE.UGZX.UGZX.DAY8 4.95 
C17 .TMKH.UGZX.QFQE.UGZX.UGZX.DAY13 8.00 
C17 .TMKH.UGZX.QFQE.UGZX.UGZX.DAY17 8.00 
C17 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY11 8.00 
C17 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY14 8.00 
C17 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY15 8.00 
C17 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY16 48.00 
C17 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY18 8.00 
C17 .TMKH.UGZX.TYFR.UGZX.UGZX.DAY19 8.00 
C17 .TMKH.UGZX.WWYK.UGZX.UGZX.DAY12 8.00 
C17 .TMKH.UGZX.WWYK.UGZX.UGZX.DAY21 16.00 
C17 .TMKH.PKVV.QFQE.UGZX.PKVV.DAYS5S 2.00 
C17 .TMKH.PKVV.QFQE.UGZX.PKVV.DAY8 35:05 
C17 .TMKH.PRVV.TYFR.UÜGZX.PEVV.DAYS 6.00 
C141.TMKH.FFTJ.TYFR.EFIJ-EPFIJODAYT23827:38 
C141.TMKH.FFTJ.TYFR.FFTJ.FFTJ.DAY27 54.42 
C141.TMKH.FFTJ.WWYK.FFTJ.FFTJ.DAY20 52.82 
C141.TMKH.FFTJ.WWYK.FFTJ.FFTJ.DAY24 30.02 
C141.TMKH.UGZX.QFQE.UGZX.UGZX.DAY3 17.77 
C141.TMKH.UGZX.QFQE.UGZX.UGZX.DAY4 2.22 
C141.TMKH.UGZX.QFQE.UGZX.UGZX.DAY6 27.89 
C141.TMKH.UGZX.QFQE.UGZX.UGZX.DAY20 41.25 
C141.TMKH.UGZX.QFQE.UGZX.UGZX.DAY22 27.24 
C141.TMKH.UGZX.TYFR.UGZX.UGZX.DAY20 24.89 
C141.TMKH.PKVV.QFQE.UGZX.PKVV.DAYS5 17.77 
C141.TMKH.PKVV.QFQE.UGZX.PKVV.DAY26 5.58 
C130 .TMKH.FFTJ.QFQE.FFTJ.FFTJ.DAY18 22.45 
C130 .TMKH.FFTJ.TYFR.FFTJ.FFTJ.DAY24 35.91 
C130.TMKH.FFTJ.WWYK.FFTJ.FFTJ.DAY20 3.93 
C130.TMKH.UGZX.QFQE.UGZX.UGZX.DAY3 20.00 
C130.TMKH.UGZX .TYFR.UGZX.UGZX.DAY7 13.12 
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.DAY2 
.DAY3 
.DAY2 
.DAY2 
.DAY6 
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.DAY6 
.DAY6 
.DAY6 
.DAY6 
.DAY6 
. DAY4 
.DAY5 
.DAY6 
.DAY6 
DAY? 
DAY? 
DAY? 


SECONDS 


USER: Operations Research Department 


Naval Postgraduate School 


ΕΙΤΕ SUMMARY 


RESTART /home/limt/final.gO? 
INPUT /home/limt/alrep.gms 
OUTPUT /home/limt/alrep.lst 
STEP SUMMARY: 0.260 STARTUP 
0.060 COMPILATION 
31.990 EXECUTION 
0.000 CLOSEDOWN 
32.310 TOTAL SECONDS 
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